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SYNOPSIS 


POSTBUCKLING ANALYSIS OP ISOTROPIC AND ORTHOTROPIC 
RECTANGULAR PLATES SUBJECTED TO EDGE COMPRESSION/SHEAR 
USING FINITE ELEMENT METHOD 

(A thesis submitted in partial fulfilment of the requirements 
for the degree of Doctor of Philosophy by PARAS RAM SINDHI, to 
the Department of Civil Engineering, Indian Institute of Tech- 
nology, Kanpur, India). 

The postbuckling analysis of a plate requires the solution 
of a set of non-linear partial differential equations. The 
simplest form of such equations is that due to von Karman. The 
closed-form solution of these eqpiations is rather tedius, if 
not impossible. In past, attempts have been made to solve them 
using different methods, for example, doxible Fourier series 
solution, perturbation technique, dynamic relaxation method, 
energy balancing technique, Galerkin method, finite difference 
and finite element methods* In some cases it has been foimd conve- 
nient to solve von Karman' s equations in terms of three displace- 
ment functions while many investigators have preferred to solve 
these equations in their modified form involving only two vari- 
ables, i-e. the stress function and the lateral deflection of 
the plate. 

In this thesis, the postbuckling analysis of simply supported, 
rectangular, thin, flat plates is investigated in considerable 
detail. Two types of load are considered: longitudinal compre- 
ssion and edge shear. The basis of the formulation is limited 
to linear elastic plates only. Both, isotropic as well as ortho- 
tropic plates are analysed using finite element method* The 
finite element model used is a 4-noded rectangle with 6 degrees 





2 

of freedom, viz. u,v,w, dw/dx, dw/dy and d w/dxdy at each node/ 
this makes a total of 24 degrees of freedom for an element. 

Only a quadrant of the plate is analysed in the case of longi- 
tudinal compression (due to symmetry vdth respect to loading 
and geometry) while the full plate has to be analysed in the 
case of shear loading* A 16-element mesh is used to discretize 
the quadrant of the plate or the plate itself, as the case may be. 
The finite element matrix equations are derived by minimizing the 
total potential energy that takes into account the stretching of 
the middle surface of the plate as well/ in the expression of 
strain energy, therefore, non-linear strain-displacement relations 
are used# Stiffness coefficients for linear and non-linear stiff- 
ness matrices are obtained as second partial derivative (with 
respect to nodal degrees of freedom), of the corresponding strain 
energy# The incremental method, supplemented with the Newton- 
Raphson iterative technique, is used to solve the resulting set 
of non-linear algebraic simultaneous equations. 

For plates xinder longitudinal compression, two types of 
inplane boundary condition on vinloaded edges are examined : 

(i) completely free (to move in its own plane), (ii) completely 
restrained (i.e* not allowed to move in the transverse direction) . 
Plates with different aspect ratios and buckling in different 
modes are studied in detail. To study the effect of orthotropy, 
two materials, viz. the glass-epoxy and boron-epoxy are conside- 
red. Results such as critical load curves, load-deflection 
curves, load-shortening curves, changes in buckle pattern, 
effective width curves, and stress distributions are reported. 



It is shown that the plate behaves quite differently under the 
two inplane boxindary conditions* 

A simple method is suggested to predict whether or not a 
changeover from one buckling mode to another takes place. 

Simple formulae are also devised to ascertain the magnitude of 
the possible secondary buckling load, under each inplane boundary 
condition. The effective-width curves presented are unique in 
that they account for the change in the buckle pattern of the 
plate, and are thus valid even after the plate has changed over 
to the higher mode of buckling. 

To study the postbuckling analysis of a plate subjected 
to pure shear, the NASTRAN model is modified and used in the 
present work. Plates of different aspect ratios are analysed 
and resxilts such as critical shear loads, load-deflection curves, 
deflection contoxirs, load-shear displacement curves and stress 
distribution are reported. It is shown that a shear loaded 
plate buckles in a mode which is essentially a combination of 
either some uncoupled symmetric modes or some uncoupled anti- 
symmetric modes, and that the changeover from a symmetric to an 
anti -symmetric mode (or vice versa) is not possible- The reduc- 
tion in the shear stiffness (from the prebuckling to the post- 
buckling stage) is very low as compared to the reduction in 
the axial stiffness for a longitudinally compressed plate. It 
is foiind that the tensile membrane stresses (longitudinal as 
well as transverse) get developed which tend to stiffen the 
plate. 


xxii 


A brief, yet complete, review of the existing literature 
on the subject of postbuckling has been presented. Wherever 
possible, the results of the present investigation are compared 
with those already available. 



CHAPTER 1 


INTRODUCTION 

Plates are extensively used in aircraft construction, 
shipbxiilding, plate girders and box girders used in buildings 
and bridges, and in many thin-walled structures. In recent 
years, the availability of high strength alloys and composites 
and the use of optimization techniques in designs for the minimum 
weight, have resulted in structures that are inherently very thin^ 
This, in turn, has enhanced the possibility of buckling. 

1.1 POSTBUCKLING - A GEOMETRICALLY NON-LINEAR PROBLEM 

1.1 •! Buckling and Postbuckling Phenomena 

A plate may be subjected to any type of inplane loads and i 
load-combinations. Under these loads, a thin flat plate initially 
undergoes inplane displacements only- With increasing loads a 
stage is reached when the flat equilibrium configuration becomes 
unstable and the plate passes to a buckled configuration, 
(characterized by the out of plane deflection). The load at 

i 

which this bifurcation of equilibrium takes place is called the | 
critical load (or the buckling load). 

. '! 

Unlike a column for which the critical load is, more or i 

i 

less, the ultimate load, the plate does not fail immediately at 
the critical load. It is well established that plates have 
considerable reserved strength even beyond the classical buckling 



i.e. in the postbuckling range. From Figure 1.1 it is noted that 
at the onset of buckling, the plate has zero flexural stiffness 
but it increases with the increase in lateral deflection. Thus, 
the plate is able to resist greater loads. The nature of the 
load-deflection curve for the plate in the postbuckling range is 
of 'stiffening' type - stiffness increasing with the load (as 
against the' softening' type for shells in which the stiffness 
decreases with load) . 

In general, plates are not perfectly flat but have sane 
initial geometric imperfections. The dashed curve in Figure 1.1 
is for an imperfect plate. The out of plane deflection grows 
continuously with the load, right from the beginning; the curve 
coincides with that of a perfect plate only at higher loads. 

The determination of the critical load for a particular 
structure is essentially an eigen-value problem based on linear 
(small deflection) theory. This yields only the critical load 
and provides no information about the magnitude of the lateral 
deflection at the onset of buckling and beyond. Therefore, the 
use of non-linear equations, based on the large deflection theory 
becomes essential for studying the behaviours of the structure 
in the postbuckling regime. 

1.1.2 Solving a Geometrically Non-Linear Problem 

The formulation of a geometrically non-linear problem, in 
general, is based on the following two approaches : 

(i) Total La gran gi an Approach 
(ii) Updated La gran gian Approach (or Eulerian Approach). 
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In total Lagrangian approach the formulation is always 
based on the initial undeformed configuration whereas in the 
updated Lagrangian approach, the problem is formulated with 
respect to the current deformed configuration of the structure# 
The total Lagrangian method is advantageous since the initial 
configuration of the structure remains fixed, thus resulting 
in the simplification of the formulation and the ensuing 
computations • 

The finite element formulation of a non-linear problem is 
usually done in one of the following ways : 

(i) The governing non-linear differential equations are 

converted into a matrix form, using Galerkin's weighted 
residual approach# 

(ii) The equations of equilibrixom are obtained by using the 

principle of stationary potential energy i#e. by setting 
the first variation of the potential energy with respect 
to the nodal displacotients, to zero. 

In both the cases the resulting non-linear simultaneous 
algebraic equations are solved by some iterative schemes, like 
Newton-Raphson method, as discussed in Chapter 2. 

1 .2 EFFECTIVE WIDTH CONCEPT 

% 

When a plate is loaded in xmiaxial compression, membrane 
stresses are developed within it# So long as the load P is less 
than the critical load P„ , the longitudinal stress, i.e. a , 


is uniform at all points, (Figure 1.2a)# Obviously at this stage. 
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Fig. 1.2 Effective width concept 





b 


P = b t = b t = b t (X.l) 

where a is the average lonaitudinal stress, o the maximTim 

stress, b the width and t the thickness of the plate. The 

stress distribution continues to be uniform until P = P when 

cr 

djj approaches (the critical stress) • For loads beyond 
tensile membrane stresses (in the transverse direction) are 
developed within the plate, which retard the increase of the 
lateral deflection. The stresses on loaded edges, therefore, 
get redistributed in a manner that more and more stress is 
produced in the vicinity of the unloaded edges and the stress 
in the central portion remains low ~ generally less than f 
(Figure 1.2b). The nonuni forniity in the stress pattern is 
actually due to the differential plate shortening along the 
longitudinal direction. 

At any stage of loading (P > P^^) the total compressive 
load may be assumed to be carried by the two strips of equal 
width, adjacent to the unloaded edges (Figure 1.2c), the central 
portion remaining ineffective in resisting the load. The 
effective width b^ is, thus, defined as 


P = t = a_„ b t 

max e av 


( 1 . 2 ) 


i .e. 


b 

/ dy 


a ^ ■ "x 

max o 


(1.3) 


The concept of effective width is useful for design 
purposes; the ultimate strength of the plate is readily estimated 


by requiring not to cross the yield stress of the material. 



From time to time, many empirical and semi-empirical 
effective width formulae have been proposed, the earliest being 
that of von Karman • A reference to these formulae can be made, 
for example, in Lind et. al. (1976) and Narayanan and Chow (1935) 

1.3 BOUNDARY CONDITIONS 

In this section the standard boundary conditions are 
summarized for rectangular plates in terms of flexural (out of 
plane) displacement/force and membrane (inplane) displacement/ 
force. 


1.3.1 Flexural Boundary Conditions 


( a ) Simply supported edges 

The deflection and the bending moment along these edges 
is zero, i.e. 


(b) 


(c) 


w = 0, M = 0 along x = constant 

(1-4) 

w = O, My = O along y = constant 
Clamped edges 

Along such edges the deflection and the slope is 
required to be zero, i.e. 

w = 0, 3w/3x = 0 along x = constant 

(1.5) 

w = 0, 9w/3y = 0 along y = constant 
Free edges 

On free edges the shear force and the bending moment 


are zero, i .e 



0, ~ along x = constant 


( 1 . 6 ) 


Qy = 0, My = 0 along y = constant 

1.3.2 Membrane Boundary Conditions 

( 1 ) Specified displacements 

The inplane displacement on an edge may be specified 
as 

u = u* along x = constant 

(1.7) 

V =s V* along y = constant 
( i i ) Specified normal membrane force 

Along an edge, some normal manbrane force of known 
value may be specified. 

( 1 i i ) Shear-free edges 

For an edge to remain shear free, the shear stress/ 
force can the edge is assumed to be zero, such that 

9 v/ 3 x + au/ay = 0 along the edge (l.S) 

(iv) Stress-free edges 

For an edge to remain stress free, both normal (direct) 
and shear stress/force on it are assumed to be zero. 

The inplane boundary conditions for the straight edge (SE), 
movable edge (ME) -which is similar to the stress-free edge and 
fixed edge (FE) are separately dealt in Chapter 3, section 

3.1. 



1 .4 


LITERATURE REVIEW 


1 *4 •! Isotropic Plates Under Edge Compression 

Levy (l942) and Coan (1951) are amongst the first few who 
initiated the study of the postbuckling behavior of rectangular 
plates by presenting approximate solutions of von Karman's large 
deflection equations* Yamaki (i960) extended the work of Levy 
and Coan by incorporating different inplane and flexural boundary 
conditions* A similar analysis for the case v/hen unloaded edges 
of the plate are rigidly held apart, has been reported in Timoshenkc 
and Gere (l96l) and has also been independently attempted by 
Rushton (1969). Brown and Harvey (1969) studied the behavior of 
plates under the combined action of lateral and edge loads* The 
perturbation technique to investigate postbuckling problems was 
employed by many investigators, for example, Stein (1960), VJalker 
(1969), Pbodes and Harvey (1971 , 1977), Dombourian et* al* (1976). 
Among other methods, Rushton (1969, 1970) used the dynamic 
relaxation technique and Brown. and Harvey (1969) applied the finite: 
difference procedure to solve the goveiming differential equations ! 
of the related problem. ; 

The Finite Element Method (FEM) to solve non-linear plate 
problems has also been employed by many investigators/ most, 
notable among then are Kawai and Yoshimura (l969), Murray and 
Wilson (1969), Yang (1970, 1971), Roberts and Ashwell (1971), 
Colville et* al* (1973), Shye and Colville (1979), Yang and Han 
( 1983) , Camoy and Hughes (1983) and Szilard (1985)* While sane 


XU 


of these investigators (Kawai and Yoshimura, 1969/ Colville et« 
al., 1973 and Shye and Colville, 1979) have used the direct 
iteration procedure to solve the resulting finite element matrix 
equations, there are some (Yang, 1970, 1971/ Yang and Han, 19 S3) 
who used the inciremental approach* Roberts and Ashwell (1971) 
made use of the mid-increment stiffness matrix in the incremental 
method and Szilard applied the energy balancing technique to the 
problem. The majority of researchers have used the total Lagrangian 
approach in their formulations, except a few like Murray and Wilson 
(1969) who adopted the updated Lagrangian approach. The finite 
element models employed in these analyses vary from simple triangles 
(Murray and Wilson, 1969) and simple rectangles (Yang, 1970) to 
Isoparametric shell elements (Carnoy and Hughes, 1983) and high- 
order triangular elements (Yang and Han, 1983) . In studying the 
postbuckling behavior of isotropic plates, some authors have also 
examined the problem of secondary buckling and of change in the 
wave-form; these include Stein (i960), Chang and Masur (1965), 

1964 . ^968 

Rushton (1970), Supple yl970) , Shye and Colville (1979), Nakamura 
and Uetani (1979) and Camoy and Hughes (1983). 

1 .4 .2 Orthotropic Plates Under Edge Compression 

A very early attempt in this regard is due to Wittrick (1952), 
who brought out a correlation between some stability problems for 
orthotropic and isotropic plates under uniaxial and biaxial 
compression. Basu and Chapman (1966) analysed plates subjected to 
uniformly distributed lateral load whereas Aalami and Chapman (1969) 



included the effect of edge loads as well* Chia (1972) and 
Niyogi (l973) used the perturbation and the Galerkin techniques 
respectively, to solve governing equations for the case of 
uniformly distributed load* 

The available work which deals with the postbuckling 
behavior of orthotropic plates is rather small « Chandra and 
Raju (1973) and Prabhakara and Chia (1973) used the perturbation 
method and the double series solution respectively, to Investigate 
this problem. In their investigations, unloaded edges of the plate 
were maintained straight throughout* Harris (1975) investigated 
both the buckling and the postbuckling of a class of biaxially 
loaded laminated plates in the presence of bending-extensional 
elastic couplings* Stein (1933), in an extension of his earlier 
work on isotropic plates (Stein, i960) has analysed simply supported 
and clamped long rectangular orthotropic plates to study the 
postbuckling phenomenon* 

1 *4 *3 Plates Under Shear Loading 

A review of the work on shear buckling of isotropic and 
orthotropic plates has been presented by Johns (1970) in a 
comprehensive style* The problem of postbuckling analysis was 
first undertaken by Kaminski and Ashton (1971) who presented 
an experimental study of the post shear-buckling diagonal tension 
behavior of rectangular laminated boron— epoxy plates clamped on 
each edge. Kobayashi et- al. (1981) also made an experimental 
study on graphite-epoxy laminated plates; the ultimate load was 



found to be considerably greater than the initial buckling load. 
Agarwal (1981 ) examined the postbuckling behavior of multibay 
composite shear webs to study failure modes typical of composite 
panels. The investigation of Zhang and Matthews (1934) on the 
postbuckling analysis of anisotropic plates under combined 
compressive and shear loading revealed the importance of the 
direction of the applied shear force in the case of anisotropic 
plates. Very recently, Stein (1985) has analysed long orthotropic 
plates in ccxnbined shear and compression; the method used is 
similar to the one developed in his earlier investigations (1960, 
1933) dealing with longitudinally compressed plates* 

1.5 OBJECT AND SCOPE OF PRESENT INVESTIGATION 

From the foregoing literature review, it is clear that 
much of the work on postbuckling has been for longitudinally 
compressed isotropic plates* In most of these investigations, 
an important feature, namely the "change in buckle pattern" 
beyond the primary buckling, has been overlooked. The resulting 
curves for the effective width and the maximum longitudinal stress 
are, therefore, applicable only for the primary mode. All the 
analyses for orthotropic plates also suffer frcxn the similar 
drawback. Regarding the consideration of inplane boundary 
conditions, the unloaded edges of the plates have mostly been 
assumed to be held straight (SE condition) . Some investigators 
have used stress-free unloaded edges that are free to move in the 
plane of the plate (ME condition) and a very few have considered 
the edges rigidly held apart (FE condition). Finally, the 


postbuckling behavior of plates under sh^ar loading has received 
very little attention; the most comprehensive work is due to 
Stein (1935) which, however, deals with long plates only# 

The literature review also reveals that the finite element 
method has been employed only for isotropic plates under edge 
compression. The method has# so far, not been extended to include 
orthotropic plates and the plates loaded in shear# With these 
points in view, the aim of the present investigation is - 

i) to demonstrate the potentiality of the finite element 

method to handle the problem of postbuckling of plates, 
taking into consideration the variation in material 
properties (isotropic and orthotropic) , load conditions 
(compression and shear) and inplane (membrane) boundary 
conditions (ME and FE conditions) • 

ii) to explore in detail the behavior of a plate xonder FE 

condition (which relatively has received less attention 
in the literature) and to compare the results with those 
under ME condition. 

iii) to study the phenomena of the change in buckle pattern, 
beyond the primary buckling. 

iv) to arrive at some formula /method to estiraate/predict 
the load at which a changeover from the primary mode 
may take place. 


v) to present the effective width curves which take due 
consideration of the change in the buckling mode. 



vi ) to get an insight into the postbuckling behavior of 
isotropic and orthotropic plates subjected to pure 
shear* 

Thus in brief, the thesis deals with the postbuckling 
analysis of simply supported thin isotropic and orthotropic 
rectangular plates subjected to either uniaxial compression or 
to pure shear. The entire thesis is divided into six main 
chapters. Chapter 1 is the Introduction. In Chapter 2, 
complete details of the finite element formulation are given, 
the governing equations are derived and the solution procedure 
is discussed. Chapter 3 deals with the analysis of isotropic 
rectangular plates subjected to longitudinal compression. This 
analysis has been extended in Chapter 4 to Include orthotropic 
materials,* the results presented are for glass-epoxy and boron- 
epoxy plates# Chapter 5 is devoted to study the postbuckling 
behavior of both isotropic and orthotropic plates under pure 
shear. Although, the results are presented and discussed at 
the end of each chapter, the main results and conclusions have 
been summarised in Chapter 6 for the sake of completeness. 



CHAPTER 2 


FINITE ELEMENT FORMULATION AND SOLUTION PROCEDURE 


2 .1 TOTAL POTENTIAL OF THE FINITE ELEMENT 


Consider a rectangular finite element shown in Figure 
2»1. The total potential energy Ji of such an element consists 
of its strain energy U and the work done V by various forces 
acting on it, so that 

jr = U - V C2 .1) 

2.1.1 Strain Energy of an Isotropic Element 


The strain energy U for a thin plate is given by 
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(2 . 2 ) 


where a , o , x are the stresses and e y the strains. 

The stresses associated with the z-direction, being very small, 
have been neglected. 

For an isotropic material, the stress-strain relations 

are 
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(2 .3) 






where E is the modulus of elasticity and v the Poisson's ratio 
Substitution of Eqs* (2*3) in Eq« (2 •2) yields 
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Et 


+ 2 6^ + kl-v) r2 ] dx dy (2.4) 
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where t is the plate thickness. At any point, the strain- 
displacement relations from large deflection theory are. 
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xy ay 3 x ax- ay- -- gxay 

in which u,v,w are the displacements in x-,y-, z-directions 
respectively. 

Substituting Eqs . (2.5) into Eq. (2*4), the strain energy 
expression becomes 


U 

where. 


+ U, + U-, 
0x2 


u „ + u ,+u, +u~ 

o,ra o,b 1 2 


(2 . 6 ) 


U 


U V. 
o,b 


- ^ < 1 ? - ^ 

D ;/ ((a!w)2 2^ (1!«){2!|) + (1^)2 

3x 3X ay 3y 

* 2(1-’') C^)"! dx <jy 


(2 .7) 


(2 . 8 ) 





37 


^2 = 


3D 


// {(|H)2 ^. (^)2j2 ^y 


(2 .9) 
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and D is the flexural rigidity of the isotropic plate, given by 


D = 


Et' 


12(1-^^) 


(2 . 11 ) 


In the above expressions, represents strain energy 

due to inplane stretching (membrane) and ^ due to bending, 
both being of quadratic order. The part represents the cubic- 
order strain energy due to coupling between the stretching and 
the bending effects- The quantity U 2 is of quartic order repre- 
senting the effect of rotations due to large deflections. 


2 .1 .2 Strain Ehergy of an Orthotropic Element 


It is assximed that this ortho tropic element has its 
material axes of symmetry parallel to the coordinate axes and 
thus the following derivations apply to a specially orthotropic 
element- The expression for strain energy U, Eq. (2-2) is valid 
for an orthotropic element also. But the stress-strain relations 
for this case are given by 
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Here and are, respectively, moduli! of elasticity 

and Poisson ratios associated with x- and y-axes and G is the 

xy 

shear modulus associated with x- and y-axes* 

The substitution of strain-displacement Eqs. (2 »5) in 
Eqs • (2*12) and then in Eq. (2*2), results in the following 
expressions for various components of strain energy for the 
ortho tropic element : 
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Here, flexural rigiditives D ,D and shear rigidity D are given 
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2 .1 «3 Work Done by Forces Acting on the Element 

Work done by inplane forces My and is 

V, = <t)u(N n + N n )dB 

1 ^ XX xy y 


(2 -18) 


V 2 = ^v(N^y n^ + Ny n^) dB 

where n^ and Hy are the direction cosiKes of the outer normal 


and the line integral is along the boundary of the element- The 
forces N^,Ny and N^y in terms of displacements for isotropic and 
orthotropic plates, Chia (1980), are : 
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The work done by lateral forces Q ,Q and moments 
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where the forces and moments, for isotropic and orthotropic 
cases (Chia, 1980) are : 
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Finally the work done by the surface load q is 


Vg = // qw dx dy (2 .24) 

2.1.4 Total Potential Energy 

Collecting various expressions for strain energy and 
the work done by loads and then substituting them in Eq. (2.1), 
the total potential energy jr for an isotropic element becomes 


71 = 


. 6D 


// 




dx dy 



SI +2(l->^l(||^)^}dx dy 

ax a? ay^ a? ^=‘2^ . 


+ 



JJ 




*^ss dx dy - q // w dx dy 

+ V V 

- ‘^''‘“xy ”x * «y "y> ® 

- SwCO^ "x ■'■ ^y °y* 

- ^a^ ‘«x "x * «xy V 

- f>ly ‘V “x-^ “y V “ 


(2 .2 5) 



The corresponding expression, for the orthotropic element is 
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2.2 DISPLACEMENT FUNCTIONS 

To guarantee mono tonic convergence, the displacement 
functions chosen should satisfy the compatability and complete- 
ness requirements. (Huebner and Thornton, 1932) . For compati- 
bility, it is observed from Eq. (2.25), [ Eq, (2.26) for ortho- 
tropic case 1, that the highest order of derivatives of displa- 
cements u,v and w, under the double integral sign are one, one 


and two, respectively# Thus compatibilities of u,v,w, dw/dx 
and dw/^y are needed# For completeness, the complete Eq.(2»25), 
[sq. (2 #26) for orthotropic case] is viewed# It is found that 
the highest order of derivatives of u, v and w are one, one and 
three respectively. Hence, for inplane displacements u and v, 
the completeness upto first order derivatives, and for the out 
of plane displacement w, the completeness upto third order deriva- 
tives are required# 

The requirements of compatibility (of zeroth order 
derivatives) and completeness (of first order derivatives) of 
displacements u and v are identical# Thus the approximating 
polynomial displacements for these will also be identical. The 
bilinear displacement rectangular element i.e. the four noded 
rectangular element with u and v as the degrees of freedom will 
satisfy the above requirements/ thus displacement functions are 

u^®^ = N2 N3 N^J 

(2 #27) 

v ^®^ = ^2 ^3 ^ 4 - 

where the superscript (e) denotes 'over an element' and (ne) 
denotes 'at the nodes of the elenent' # Interpolating functions 
to can be found in Zienkiewicz (1977) # 

The requirements of compatibility and completeness of 
displacement w are more stringent than those of u and v# Compati- 
bility upto first order derivatives and completeness upto third 
order derivatives are needed. The bicubic rectangular element 



i.e. the fibtir noded rectangular element with w, ~ , ~ and 
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degrees of freedom at each node will satisfy the above 
requirement. Thus displacement function w is 
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Interpolating functions to can be found in Zienkieicz 
(1977). Thus the rectangular finite element model. Figure 2.1, 


will have a total of 24 degrees of freedom with 6 at each node 
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2.3 BASIC EQUILIBRIUM EQUATIONS 

2.3*1 Equilibriiim Equations For the Element 

Using Eqs. (2*27) and (2*28) the expression for total 
potential energy [Eg. (2-25) or Eq. (2*26)] , becomes 
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where, now, U^, Uj^ and U 2 are, respectively, quadratic, cubic 
and quartic functions of nodal displacements 
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The equations of equilibrium are obtained by minimising the 
above energy ; 
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where is the load matrix of external loads and {B} ^*^®^ 


is the boundary force matrix 



Equations (2 *30) can be put in the form 
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It may be noted that the stiffness matrices 
and [n 21^®^ are all symmetric. The matrix is the 

conventional linear stiffness matrix and [ ] and [N 2 ] are non- 

linear stiffness matrices which are respectively linear and 
quadratic functions of nodal displacements. Many a times Eq. 
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is also 

written as 
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where 
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The matrix is called the secant stiffness matrix . 

2 . 3.2 Evaluation of Matrices 
Case I t Isotropic Plate 

For convenience in evaluating these matrices^ the vector 
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{ q} ^ ■' is rearranged as 
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With this in view, tte matrix can be expressed as 

Various submatrices of Eq. (2.36) are given by 
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In the derivation of above expressions# Eqs . (2*7), (2*8), (2*27) 
and (2*28) have been used* 


The non-linear matrix 



is similarly expressed 


as 



where 


l^"l,23l 1^^1,331 


n 


1,13 


3{u} 


(ne) 


3{w} 




(2 .41) 


^ // [2£N,x} Lm,x J LM,xJ 


+ 2V {N,x} LM^yJ LM,yJ 

+ 2(1-V) £N,y} LM,xJ L M,yJ ] dx dy (2*42) 

1,2 3 3{v}^”®^ 3£w} ^’^®^ 

= —2 // [2^^ £N,y} L M, xJ £w}^’''®^Lm, xJ 
t^ 

+ 2£N,y}LM,yJ £ w} ^’^®^ LM,yJ 


+ 2(1-^) £N,x} L M,yJ £w} ^”®^ LM,x J ] dx dy (2*43) 


n 


1,33 


3£w} 


a^u. 

(neT 


a£w} 


rneT 



LM,x_I 


= “ j / / [2{M,x} L N#xJ {u} 

+ 2V {M,x} LN,yJ {v}^"®^tM,xJ 
+ 2{M,y} LN,y J {v} t_M,y J 

+ 2V {M,y} LN,xJ{u}^"®^ LM,yJ 

+ (l-V){M,x} Ln,xJ Cv}^”®\LM,yJ 
+ (l-a^){M,y} LN,xJ {v}^”®^ LM,x J 
+ (l-5^){M,x} LN,yJ {u}^''®^ LM,yJ 

+ (l“V)£M,y} LN,yJ {u}^’^®^ LM,xJ ] dx dy (2-44) 


In the derivation of above relations, Eqs . 
(2 #28) have been used. Finally, the second 
is expressed as 


(2 #9), (2 #2 7) and 
non-linear matrix 




where 


0 



0 O 

° K.33I 


(2 #45) 


”2,33 


= “1^ [{M,x} |_M,xJ Cw}^”®^ LM,xJ {w}^”®^ LM,xJ 

+ {M,y} LM,yJ {w} ^"®^ LM,yJ { w} ^’^®^ LM,yJ 



+• {M,X} LM,xJ {w} LM/Y-J {w} LM,y J 

+ {M,y} LM,xJ LM,yJ L M, x J ] dx dy 

(2.46) 

Equations (2.10) and (2.28) have been used to derive the above 
matrix equation . 

Case II ; Orthotropic Plate 

Submatrices of and [n 2 ]^®^ for the 

orthotropic case are similarly derived using the corresponding 
expressions of strain energy {Eqs. (2.13) - (2*16)./] and are 

SJ (2D^{N,x} LN,xJ+ 2 D^{N,y} LN,yJ.)dx dy (2.47) 
"t 

k 3^2 = SS (2D^Vy{N,x} LN,yJ + 2 D^{N,y} LN,xJ)dx dy 

(2 .48) 

^22 “ "T (2Dy{N,y} LN,yJ + 2 D^{N,x} LN,xJ)dxdy (2.49) 

k33 = // (d^{M,xx} LM,xxJ+ Dy{M,yY} LM,yyJ 

+ D^Vy{M,xx} LM,yyJ+ Dy {M,yyl LM,xxJ 

4 LM[,xyJ)dx dy (2.50) 

xy 

n^ ^3 * J/ (2 D^{]:J,x} Lm,xJ £w}^*^®^ LM,xJ 
^ "t 

+ 2 D^-Vy{N,x} LM,yJ L.M,y-J 

+ 4 D fN,Y} LM,x J {w} ^”®^ LM,yJ) dx dy (2.5l) 





”l,23 = (2 DyV^{N,y} LM,xJ {w} LM,xJ 

+ 2 Dy {N,y} LM,yJ {w} L M,yJ 
+ 4 D^ {N,x} LM,yJ Lm,xJ) dx dy (2.52) 

^1,33 D^^CM^x} L N,xJ LM,xJ 

+ 2 D^Vy {M,x} LN,yJ LM,xJ 

+ 2 Dy {M,y} LN,yJ {v}^''®^ LM,yJ 
+ 2 DyV^ £M,y} L N^xJru}^"^®^ LM,yJ 
2 CM,x} LN,xJ {v}^^®^ LM,yJ 

+ 2 D^ £M,y} tN,x J£vJ^”®^ LM,x J 
+ 2 Dxy £M,x} LN,yj£u}^”®^LM,yJ 

+ 2 £M,y} L N,y J £u}^^®^ L M, x J ) dx dy (2 .53) 

”2.33 = ^ ^°x Lm.xJ £w}^”®^ Lm.xJ£w}^”®^ LM.xJ 

+ Dy £M,y} LM.yJ £w}^”®^ LM.yJ £w}^”®^ LM,yJ 
+ 2 D^Vy £M,x} LM,xJ£w}^”®^ LM,yJ£w}^”®^ LM,yJ 
+ 4 £M,y} LM.xJ £w}^”®VLM,y J £w} ^”®^ LM.xJ ) dx dy 


(2 .54) 


2 *3 *3 Equilibrium Equations for the Plate 

Elemental equations (2 •31) are assembled in the usual 
way, and on application of boundary conditions, take the foll- 
owing form 

^[^ol I + I [N2I) (2.55) 

or [Kg] (2 .56) 

and. are simply the simultaneous non-linear, algebraic equations. 
Here the superscript (n) denotes 'at all the nodes of the plate',* 
n thus represents total number of nodal displaconents . 

2 .4 SOLUTION PROCEDURES 

The equilibrium equation (2.56) is a set of non-linear, 
simultaneous, algebraic equations for the complete plate, which 
must be solved. Mainly, three methods, viz. the direct iterat- 
ion, the Newton-Raphson and the incr«nental method are used to 
solve these equations (Zienkiewicz, 1977; Cook, 1981; Yang, 1970; 
Stricklin et. al, 1971). 

In the direct iteration method the full load is applied 
in one step. To start with, some solution is guessed (typically 
{q}^ = {0}), and Eq. (2.56) is used to find { q) . Then {<!}]_ is 
used to update [Kg], and {q }2 is determined, and so on till the 
convergence is obtained. The main drawback of this method is 
that for the problons involving high non-linearity, the solution 
converges very slowly or may even diverge. Morever, the method 
is computationally expensive. 



In the Newton -Raphson method when an approximate 
solution { for a given load { Q} is reached, Eq. (2.56) is 
not satisfied, 

(2 .57) 

Let the improved solution { + Aq} makes the above inequality 
to vanish. 


^ (2.58) 

Using first order Taylor expansion, this yields 

[Kg] {q^} + ([Kg]{q})lg {Aq} - { Q} = {0} 

or [k^] {Aq} = {Q} - [Kg] {q^} (2.59) 

where ^[Kg] {q})lg = ([k^^] + [n^^] + [^2^^^q 

(2.60) 

and is called tangent stiffness matrix. The right hand side of 
Eq. (2.59) may be interpreted as the unbalanced force vector, say 
{ g} . This equation is solved for {Aq} and then using {q^=qQ+Aq}, 
the iterations are carried until {Aq} is sufficiently small. 

This method will converge slowly if the initial guess is 
not close to the exact solution. The method can also diverge for 
highly non-linear problems. 

In the incrgnental method , it is assumed that the exact 
solution {q^} of Eq. (2*56) is known for load {Q^}. and a solution 
{ q^ + Aq} is now desired for { Qq + AQ} ; Eqn. (2.56) beccmes 


(2 .61) 


+ Aq) {q^ + Aq} = {Qq + 

Using a first order Taylor's expansion, one gets 

[K 3 I {q}j ([Kg] [q}j {Aq} = { Q_} t {AO} 

% % 

or [k^] {Aq} = {aQ} (2.62) 

where [k^] is the tangent stiffness matrix given by Eq. (2.60). 
Thus in this method the load is applied in increments and 
corresponding displacements are found using the latest value 
of the tangent stiffness matrix. It is claimed ( Zienkiewicz, 

1977) that by using suitably small increments of load the method 
guarantees to converge. However, the method has a disadvantage 
that no real estimate of solution accuracy is known since, in 
general, the equilibrium is not satisfied at any given load level. 

To overcome the above difficulty the incremental method 
is usually supplemented with the Newton-Raphson method; in this 
combined form, it is known as the lncremental-cxim~lterative 
method and the same has been used in the present work. The load 
is applied in small steps of size £AQ} and the corresponding 
incremental displacements £Aq} is obtained from Eq. (2-62). After 
a few steps, the updated (approximate) displacements {q^} and 
load { Q} are used in the Newton-Raphson equation. Then the 
Eq. (2.59) is solved for { Aq} , and iterations are carried till 
{Aq} is reduced to zero or to a prescribed tolerance. In this 
way, almost an exact solution is obtained. The process is 
repeated till the full load is reached. 


Stress Calculation s 


Membrane strains and stresses (at the middle surface of 
the plate) are calculated from Eqs- (2-5) and (2.'3). It is well 
known that the stresses calculated at the nodes are usually in 
error, ( Zienkiewicz, 1977, Cheung and Yeo, 1979) . Further, it 
is also known that for linear elements, stresses are accura- 
tely determined at the centroids of the elements. Similarly for 
elements, best results are obtained at the Gauss points. For 
non— linear strain— displacement relations (Eqs. 2.5), the linear 
parts, therefore, need to be evaluated at centres and the non- 
linear parts at the Gauss points for obtaining accurate results. 
However, it is foiand that if both the terms (first two terms of 
Eqs. 2.5; the third term being zero at the middle surface) are 
evaluated at the centre of the element, the stresses are accurate 
enough. Hence, the strains (and stresses) are found at the 
centres of the elements and then extrapolated to the nodal points 
by bilinear extrapolation. 

Some Computational Aspects 

For the solution procedure adopted, in the present work, 
the [k^] and [Kg] matrices are obtained by numerical integration 
(Gauss quadrature) technique. Number of Gauss points used to 
evaluate [k^ ] is 3 and for [ Nj_ ] and [N 2 I it is 4. In the analysis 
[Kq] is calculated only once, whereas [n^] and [N 2 ] are updated 
for the latest displacements. 

A small lateral load is applied to initiate buckling. 

The order of this load is 0.1% to 1.0% of the critical load. 



The lateral load is maintained constant while the inplane load 
is varied (increased) to study the postbuchling behavior* 

In the prebuckling range# a load slightly less 'than the critical 
is applied in just one step. The convergence for this load is 
easily obtained in 2-3 iterations. Near the critical load# on 
either side# the size of incremental load step is kept very 
small to obtain easy convergence. Once the critical load is 
passed, the larger step size is used. The convergence is 
considered achieved when the error 

fe} = ^ ” {qii (2.6 3) 

satisfies either of the following norms for the absolute error 
and the relative error, viz. 

I e j I ^ .00001 j— l,2,...#n (2 .64 ) 

^ < .001 j = l,2,...,n (2.65) 

where n is the total number of nodal displacements. 

To accelerate the convergence in Newton-Raphson itera- 
tive method# a technique proposed by Crisfield (1979) is used. 

In this method the usual correction , corresponding to 

"til 

the current unbalanced force {g}^ , obtained in the i iteration 
(Eq. 2.59), is modified to {6}^^ then added to the current 

displacements £q}j^ to get the next (updated) displacenents 

Cqli+l 


{qli^l = iq}i + {6}i 


(2 . 66 ) 



^ / 


The accelerated displacement vector is given by 

®i ^^^i-1 (2.67) 

vrtiere is the previous accelerated displacement vector 

and and fj_ are scalar multipliers (initially e^ = O and 
f = 1 ) / given by 


fi 


where a^^ , 


- a^/b 


( 2 . 68 ) 


= (1 - Cj^/bj^) - 1 

bji^ and c^^ are scalar inner products expressed as 


b^ ({g}^ - (2.69) 

Cj^ = LAqj^ ({g}jl^ - {gj^.j^) 

in which is previous imbalanced force vector. 

2.5 TEST PROBLEM 

The computer* program developed for the present study, 
is first tested for an isotropic square plate (with side 
a s» 1 .0 m, thickness t = 10.0 mm, modulus of elasticity E = 

200 GPa and Poisson's ratio v = 0.316), subjected to uniform 
lateral pressure p, per unit of area, when the edges of the plate 
are (a) simply supported and (b) clamped. Further, the edges 
are held immovable i.e. no inplane movement normal to the edge 
is allowed. The boundary conditions, thus, are (Eqs. 1.4 and 
1.5) 
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u = w = = 0 along the pair of edges parallel 

to y-axis 

V = w = = 0 along the pair of edges parallel 

to x-axis 

for simply supported edges, and 


( 2 . 70 ) 


3w ^ 
u = w = -- = 0 


along the pair of edges parallel 
to y-axis 


(2 .71) 


V = w = = 0 along the pair of edges parallel 

to X-axis 

Due to symmetry, only one quadrant of the plate is 
analysed, which is discretized into 16 rectangular elements of 
the same size (4-noded elements with 24 degrees of freedom) . 

The consistent load matrix is used in the analysis. Load- 
deflection values are compared with those of Levy (1942), based 
on the analytical solution, and with those of Yang and Han (1933) 
who used FEM employing high-order triangular elements with 54 
degrees of freedom. The agreement amongst three results is found 
to be extremely good, as can be seen in Figure 2.2. 
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(b) Clamped edges 

Fig. 2.2 Load -deflection curves for isotropic 
square plate under uniform pressure. 





CHAPTER 3 


ISOTROPIC PLATES SUBJECTED TO LONGITUDINAL COMPRESSION 

The postbuckling behavior of isotropic plates subjected 
to longitudinal compression has been investigated by many resear- 
chers, as has been reviewed in Chapter 1* Most of these inves- 
tigations suffer from a basic drawback that the phenomenon of 
change in buckling patterns has not been incorporated* In this 
chapter, this aspect has been properly accounted for in studying 
the postbuckling behavior of isotropic plates under longitudinal 
compression . 

3*1 ANALYSIS DETAILS 

The plate considered is rectangular, of size a x b and 
with symmetric boundary conditions, and is compressed axially 
along the edges AD and BC, Figure 3»l(a)* Due to symmetry, only 
one quadrant is analysed/ the details of the quadrant APOJJ are 
given in Figure 3*l(b)* It is discretized into 16 finite elements 
of 4-noded rectangles of equal size. The plate is assumed to be 
simply supported on all the four edges • The appropriate boundary 
conditions as given by Eq. (1.4), are applied. 

Along the unloaded edges, the inplane boxindary conditions 
could be one of the following : 

(i) Movable Edges or Stress -free Edges (ME) 


u V 0 , V / O 


( 3.1) 


(ii) 


straight Edges (SE) 


This imposes a condition of constant v along edges y = 0, b 
so that 

|i=0 (3.2) 

(iii) Fixed Edges (FE) 

This implies that the edges are rigidly held apart; i.e. 
along y = 0, b one must have 

V = 0 (3.3) 

In this chapter, only two inplane boundary conditions, viz. 
the ME and the FE, have been studied. Only for the sake of compa- 
rison, one problem of a square plate is analysed by applying 
uniform displacement along x = 0, a, while the edges y = 0, b are 
kept stress-free. For all other problems uniform load has been 
applied on the loaded edges. 

Buckling Modes 

Let a particular mode be designated by (m,n), vrhere m and 
a are the nximber of half waves along the x- and y-direction, 
respectively. Since the applied load is acting only in the x- 
ii recti on, n will always be one at the critical load. 

The buckling mode could be symmetric or antisymmetric about 
:he plate centre lines . The symmetric mode requires an additional 
iondition of zero slope (9w/3x = O) across the centre line. For an 
intisymmetric mode, the corresponding condition is of zero deflec- 
-.ion (w = O) . The various modes investigated are (1,1), (2,1) and 
3.1) depending upon the problem. 


I 



To initiate a particular mode, a very small lateral load 
s applied at the appropriate position (node) . For example, for 
he mode (l,l), this load is applied at the centre of the plate 
.e. at the node 25 (Figure 3*lb) • For the mode (2,1), the 
osition of the load would be at the node 15, whereas to initiate 
(3,1) mode, three such loads are required to be applied — one 
t node 25 in some direction and the other two at nodes 10 and 15 
1 opposite directions to the previous one. It is worth pointing 
lat the disturbing force applied to initiate the desired buckling 
ittem is too small to cause any effect on the post-buckling 
ihavior. It may be maintained on the plate throughout the 
lalysis (Murray and Wilson, 1969) or removed once the load has 
•ossed its critical value or at sane other convenient point 
spending upon the method used (Colville, et, al., 1973), 


ryes for Buckling Load 

For a simply supported, rectangular isotropic plate, under 
ge compression, the critical load N (per unit width) is given 
(Chajes, 1974) 



(3,4) 


en the unloaded edges are movable (ME condition) , Here, k is 
3 compressive buckling coefficient and is given by 

’'c “ ‘■I * ' E ' 1 

i curves obtained from Eq, (3,5) are shown in Figxire 3,2* The 
insition from m to m+1 half waves occurs when 



X = '/' m(m+i ) 


(3.6) 


i.e. the successive values of X are vA 2 ... • The 

point of intersection of curves of higher modes# i-e- of m = l 
and 3 or of m = 2 and 4, is obtained from 

X = '^mTm+2T (3.7) 

the successive values being '^3 , V'8 , Vl5 ... . 

When the unloaded edges are held rigidly apart (i.e. the 
FE condition implying v = 0 along y = 0,b), the problem beccmes 
similar to the case of a biaxially loaded plate except that the 
ratio of the loads N and N is a constant and is equal to the 

y X 

Poisson's ratio ^ . The buckling coefficient k^, now is given 
by (Szilard# 1974) 

(m ^ 

A m 

k = # n = 1 (3.3) 

^ .2 2 

/ - X n \ 

(1 ^ 

m 

The curves obtained using Eq. (3*8) are shown in Figure 3.3. It 
can be seen that the transition frcm m = 1 to m = 2 and from 
m = 2 to m = 3 occurs at X = 2.48# 4*12 respectively 
(Nakamura and Uetani# 1979). 

A comparison of curves in Figures 3*2 and 3.3 leads to 
the following observations : 

i) For a particular aspect ratio# k is smaller for FE than 


for ME condition.. 


ii) The absolute minimtjin value of is 4*00 for ME (or SE 
condition) and is 2.30 for FE conditions, 
iii) A plate of fixed aspect ratio may have different primary 
modes under the two membrane boundary conditions# For 
example, for X= 2 #0, the plate buckles in (2,1) mode 
for ME conditions but in (1,1) mode under FE condition- 

3-2 ILLUSTRATIVE EXAMPLES 

Plates with aspect ratios 1.0, 1*7 and 1.9 are considered 
for the analysis# These ratios have been so selected as to lie 
in the range of (i) upto '^2, (ii) between '^2 and '/’3 and 
(iii) between '^3 and '^6, of the ciarves in Figure 3#2 for ME 
inplane boundary conditions. It is expected that this range of 
aspect ratio will cover all the salient features of the post- 
buckling behavior of rectangular plates. With this in view, 
the following examples are taken up in detail : 

I Square plate (X = l.O), 1.0m x 1 .Cm and 10 mm thick, 

E = 200 GPa, V = .316. 

(i) Unloaded edges free to move (ME condition) — (1,1) 
and (2,1) modes 

(ii) Unloaded edges rigidly held apart (FE condition) 

(1.1) arid (2,1) modes. 

II Rectangular plate (X = 1.7), 1 #7m x 1.0m and 10 mm thick, 

E = 200 GPa , V = #316 

(i) Unloaded edges free to move (ME) — (l,l), (2,1) and 

( 3.1) modes 

(ii) Unloaded edges rigidly held apart (FE) — (1,1) and 

(2.1) modes . 



Ill Rectangular plate {X= 1.9), 1 .9in x 1.0m and lOrrni thick, 

E = 200 GPa, V = ,316 

(i) Unloaded edges free to move (ME)- (l,l), (2,1) and 

(3.1) modes 

(ii) Unloaded edges rigidly held apart (FE) (1,1) and 

(2.1) modes. 

3.3 numerical results 

3.3.1 Critical Loads 

The results are reported in terms of non-dimensional 
parameters in the form of curves and tables. Wherever possible 
they are compared with the results available in literature. 

In each case of the illustrative example, the critical 
load is determined as the point of inflexion on the load- 
deflection curve or as the point of intersection between the 
pre-buckled line (the portion of the load-shortening curve, 
before buckling) and the post-buckled line (portion, after 
buckling) . Table 3.1 lists all such values in terms of buckling 
coefficients, for different cases. For the sake of comparison, 
the exact values, as obtained fron Eqs. (3.5) and (3.3) are also 
reported. As can be seen, the two set of values match beauti- 
fully. The small difference in the two values is due to the 
error in reading the inflexion point of the load-deflection curve 

3.3.2 Load-Deflection/Shortening Curves 

Figure 3.4 shows the graph between the non-dimensional load 

' ■ 2 ' 2 ■ 

N , (= N b /ti D) and the non-dimensional lateral deflection 

■ -Kl . 
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being the maximtjin lateral deflection) for a square 
plate in (1,1) mode. The curves coincide with those of Yang and 
Han, (1983) and Rush ton, (19 70) for ME condition and with Shye and 
Colville, (1979) and Uemura and Byon#(l977) for FE condition. The 
symmetric and the anti -symmetric buckled configurations, drawn on 
the longitudinal centre line of the plate, (Figure 3.5), match 
closely with Shye and Colville (1979) for FE condition. Similar 
configurations for ME conditions have been drawn in Figure 3.6. 

The - w^curves for aspect ratios 1.7 and 1.9 are presented in 
Figure 3.7* 

Figures 3.3 to 3.11, depict load-shortening curves, drawn 
between and A (= AE b t/a ji D , A being the average shortening 
of the plate at the load N^) • It is observed frcm the graphs that 
the axial stiffness of the plate has reduced to about 5096 - 70% of 
its initial value in changing from a pre-buckled state to a buckled 
one. This is the primary buckling. A further reduction in the 
stiffness takes place in passing frcra the primary mode to the 
secondary mode, when the stiffness becomes about 40% - 50% of the 
original (unbuckled) value, (Carnoy and Hughes, 1983) . Obviously, 
the magnitude of reduction is greatly influenced by the in-plane 
boxmdary conditions at the unloaded edges (the reduction in stiff- 
ness being more under ME condition than under FE condition) . The 
nature of the N - A curves is found to be in very good agreement 
with that of Rushton (l969) for FE conditions and with Rhodes and 
Harvey (1977) for ME conditions. 
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3-3»3 Change in Buckling Pattern 

There is a clear point of intersection between the primary- 
ana the secondary modes of a square plate, for both the types of 
in-plane boundary conditions (Figure 3.8). This point of inter- 
section indicates the possibility of a change in the buckling 
pattern. A similar observation was reported by Stein (i960) who 
stated that the actual changeover takes place at the load 
greater than that given by this common point. Chang and Masur, 
(1965) have differentiated these points (of intersection) from 
those at which the energy levels in the two modes become equal. 

They define the point of intersection as the point at which the 
secondary mod^ becomes stable, whereas the point at which the 
two energies become equal, gives lower bound to the secondairy 
critical load. The point of equal energies is usually at higher 
load than the point of intersection of the two N - A curves. 

For the square plate, the intersection points in Figure 
3 .8 are at = 9*70 and 7*90 for ME and FE conditions, respectively. 
Uexnura and Byon,(l977) have reported this value to be S.SO for SE 
condition/ this appears to be in order because the SE condition 
is somewhere in between the ME and FE conditions* 

A similar change in the buckling pattern is also seen in 
the rectangular plate ( X= 1*9) under ME condition (Figure 3.10b), 
and in the rectangular plates ( X = 1»7, 1*9) under FE condi-tion 
( Figure 3 .11 ) . 





3.3.4 Effective Width Curves 


Effective width is calculated using the expression 
(Chandra and Raju, 1973) 


b 


e 


N a 

X 

At E 


(3.9) 


where b is the non-dimensionalized effective width (= b /b) . 

<= e 

The curves are drawn between and b^. Curves of Figure 3.12 
are for a square plate, with both in-plane boundary conditions, 
while in Figure 3.13 the curves are for a rectangular plate 
( ^ = 1*9) . When compared with the values of b^, for the primary 
mode, obtained from the formulae used in Narayanan and Chow (1995), 
the values of present analysis (curve saob) are slightly smaller 
for the same ME condition. The same trend is found for the FE 
condition as well, when the results are compared with Timoshenko 
and Gere (l96l) . This could be due to the fact that in the present 
analysis the uniform load has been applied as against a uniform 
displacement. 


3 . 3.5 Stresses 

Stresses are calculated at the selected load levels for 
all the examples considered. Figures 3.15 and 3.16 show, respect- 
ively, the stress distribution at load N = 7.0, for a square 
plate buckling in (l,l) mode for ME and FE conditions. The 
stresses match with Coan (1951) and Rushton (1970) for ME condit- 
ions* The comparison for FE conditions can not be made due to 
non-availability of results in literature. The variation of 
0^^ (= a b^t/ir^D) with load N is plotted in Figure 3.17 for 

lUaX ni3.X X 



a square plate and in Figur^ 3 •IS for rectangular plates. 

3.4 DISCUSSIONS 

3-4 •! Symmetric and Anti -symmetric Modes 

In the postbuckling analysis by FEM, a particular mode 
is initiated in the plate by applying a very small force at an 
appropriate position (node), to have either symmetric (m = 1,3, 

5 •••) or anti-symmetric (m = 2,4/6 •••) configuration; n in 
both cases being equal to 1 • For the symmetric mode, the 
condition of zero slope (■|^ = 0) and for the an ti -symmetric mode, 
the condition of zero deflection (w = o) is set across the centre 
line of the plate. 

It is observed that if an anti -symmetric mode is initiated, 
the final converged mode is the same (initiated) antisymmetric 
mode. However, if a symmetric mode is initiated, it converges to 
a symmetric mode but not necessarily to the same initiated mode. 
Thus, a mode starting from either m = 1 or 3 ends up in a symme- 
tric mode, which could be m = 1 or 3 or its combination, at a 
particular load level. Similar observations were also made by 
Shye arid Colville (1979) . 

3 . 4.2 Changes in Buckling Modes 

An important conclusion drawn from the results is that 
one can predict in advance whether or not the change-over frcm 
one mode to another is possible. The procedure to do so is 


outlined as ; 
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1* On the graph of buckling curves (Figures 3.2 or 3.3) draw 
a vertical line through the point of given aspect ratio 
2- If it first cuts the curve of lower mode, say 'm' and then 
the next higher, i.e* *m+l'; then the changeover from 
mode m to m+1 will always take place. The corresponding 
load shortening curves would have a definite point of 
intersection . 

3. If the vertical line first meets the curve of higher mode 
and then the lower one, the changeover cannot take place. 
The corresponding load-shortening curves in such a case 
would never meet. 

Further comments have been made with regard to each inplane 
boundary condition. 

Movable Edges (ME) 

Referring to Figure 3.2, one would thus expect that 

i) for values of X upto'^2, lines of m = 1 and m = 2 would 
intersect and therefore changeover from m = 1 to m = 2 
mode is possible. 

ii) from X = '/’2 to '/’3 point of intersection of corresponding 
N - A lines is not possible. Rushton (1970) has shown 

IXl 

that for X = 1.5, the load-shortening lines of m =1 and 
m = 2 do not meet, but has not assigned any reason to it. 
This fact is further suppoirted by Figure 3.10(a), where 
for X = 1 .7, lines of m = 1 and m = 2 diverge. However, 
the changeover from m = 2 to m = 3 r^ode is possible. 
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iii) from X = ^/'3 to '/’6, again, the changeover (from m = 2 to 
m = 3) is possible because, in this range, the vertical 
line through a given X meets first the lower mode (m=2) 
curve and then the higher mode (m = 3) curve, 
iv) An interesting point was observed while executing the 

program for the plate of X = 1 . 7 . The vertical line from 
this value (Figure 3.2) cuts the various curves, such that 


crl _ ^crl * ^crl 
m=:2 m=l m=:3 


( 3 . 10 ) 


When a (2,1) mode was initiated, there was no problem in 

executing the program at any load level. However, when a 

( 3 , 1 ) mode was initiated, there was no problem in the 

beginning but as the load approached (increment by 

increment) the value P__l , the mode changed from (3,1) 

^ m=l 

to ( 1 , 1 ) and continued to ranain in it for higher loads. 

It happened so because the critical load for the first mode 
is less than that for the third mode while the zero-slope 
condition is the same for both the modes which are 
symmetrical . 

This t 3 fpe of trouble is of computational nature and is 
inherent in the incremental method. In direct iteration method 
( Shye and Colville, 1979), one may apply directly the required 
(total) load and get the solution by iterations. 

However for the plate with X = 1.9, the above problem was 
not encountered because in this case 


cr 


I < ^crl , < '■crl 


m=2 


m=3 


‘m=l 


( 3 . 11 ) 



Fixed Edges ( FE) 


A reference to Figure 3.3 reveals that 
i) for X upto 2.43, the vertical line through it meets 
different curves in a manner such that 


ii) 



< P 


cri 


m=2 



( 3.12 ) 


Hence for plates of aspect ratios upto 2 .48, the primary 
mode is (l,l) and the secondary mode is (2,1), and the 
changeover from a lower to the higher mode will certainly 
take place. Here, it is worth mentioning that Shye and 
Colville (1979), while using the same FE inplane boundary 
conditions have assumed (2,1) and (3,1) as the primary and 
the secondary modes, respectively, for a plate with X =s 2 , 
which is incorrect. 

for values of X lying between 2.48 and the following 

holds 


P 


cr 


m=2 


< 




(3.13) 


Obviously (2,1) is the primary mode in this case. But if 
(3,1) mode is initiated first, it will converge to (1,1) 
mode, as explained for ME condition. The changeover fron 
higher mode to lower is, of course, ruled out. 
for X > 4.12 the following inequality is valid 


m-3 m=x 


( 3.14) 


In this case (3,1) is the primary mode but the other two 
modes, viz. (2,1) or (1,1), being lower, the changeover is 


iii) 
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not possible. (However the changeover from (3.1) 
to ( 4 , 1 ) can take place) . 


It is thus obvious that inplane boundary conditions on 
unloaded edges play a very important role in deciding the 
possibility of the change in the buckling pattern. For example, 
as it has been mentioned earlier, for a rectangular plate with 
aspect ratio 1»7, there will be no change of mode, from (l.l) to 
^ 2 . 1 ). for ME condition but the mode does change for FE 
condition [from (l.l) to (2.1)]. 

It is further observed that the load at which the lines 
corresponding to modes (l.l) and (2.1) intersect, is greater 
for ME condition than for FE condition for the same plate : 

For example, for a square plate, the respective values of the 
load are 9.7 and 7*9 (Figures 3»8a.b). This fact is' further 
supported when the vabe 7.0. as obtained by Rushton (1970) for 
= 2 .0 plate for ME condition is ccmpared with the value 5.6 
obtained by Stein (I960) for SE condition. 

A more closer examination of results (of points of inter- 
section obtained for plates of various aspect ratios) available 
in literature reveals that for SE condition the point of inter- 
section of load-shortening cxirves. lies approximately at the 

same distance above the critical load N _ . as the later lies 

xs 

above the load N,^ . i.e. 

xp 


N.-N„ = N„ - 
XI XS xs xp 


(3.15 


where N and N are the critical loads corresponding to 

XS 

secondary and primary modes, respectively. 
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For ME condition, it was observed that the corresponding 

point of intersection is about 15-25% higher than that for SE 

condition, the higher percentage being for the larger aspect 

ratios. On the contrary, for FE condition the points lie at lower 

values than those given by the above equation • The reduction 

is small and varies between 8 - 10% only. For higher values of 

X, it is observed that the ratio of N . and N , for FE condition, 

XI xp 

is lower than that for ME condition. This was also reported by 
Nakamura and Uetani (1979). For example, for a plate with 

1.9, the ratio xinder ME condition is 1 .81 while it is 1.13 
only for- FE condition . 

3 . 4.3 Effective Width Curves 

Many empirical formulae and curves for effective width 
have been given and used in the past, they are mainly based on 
experimental work. In such formulae, it is tacitly assumed that 
the form of the buckle remains same at all stages of loading. 
However, the buckling pattern does change, (as has been discussed 
in the previous section), and it must be taken into account. 

The effective width curves of Figures 3.12 and 3.13 are, 
therefore, applicable even after the mode of the buckled plate 
has changed to the higher one. The enveloping curve saod, (with 
firm line), is to be used to estimate the effective \^ddth. Until 
buckling starts, the whole width of the plate is effective (shown 
by horizontal line sa) and thereafter, the immediate curve ao 
(corresponding to primary mode) gives b^. The effective width 
thus reduces with increase in load. At the point where the 
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second curve (corresponding to secondary rnode) cuts the initial 
curve Cat o ) , the change in mode takes place. This should be 
noted that the changeover point has (and should have) the same 
value as the point of intersection of the corresponding load- 
shortening lines discussed earlier. 

It is interesting to note that for rectangular plate 
(X = 1.9), with FE condition the changeover point is rather 
quick i.e. the mode changes from (1,1) to (2,1) immediately 
after the buckling has set in. In this case, therefore, it is 
the secondary mode which mostly dominates. This shows how 
important it is to consider the change of mode in estimating 
the effective width. 

The results of b^ versus load ratio N, for the square 
plate under ME condition, loaded with uniform displacements are 
shown in Figure 3.14(a). It can be seen that the effective 
width is greater when loaded edges are uniformly displaced than 
when they are uniformly loaded. This is true for both modes, 
i.e. (1,1) and (2,1). 

3.4.4 Stresses 

Following observations are made from the stress results 

presented in Figures 3.15 - 3.18. 

i) In the case of ME condition, the maximum longitudinal 

stress O , always occurs at the crest of the buckle 
moix 

and on the unloaded edges, whereas for the FE condition 
it develops at the junction of the loaded and the 
unloaded edge (Figures 3.15 and 3.16). 



ii) Beyond the critical load the transverse manbrane stresses 

(which were initially ccxnpressive throughout) progressively 
reduce near the central portion of the xmloaded edge and 
become tensile (Figure 3*16). 

iii) The variation (increase) of with the applied load 

is linear even in postbuckling range, but at an 
increased rate. (Figures 3*17 and 3.18). 

iv) From Figure 3.17, it is easily interpreted that for the 
same load ratio N /N^ , the stress a , is greater in 

the higher (2,1) mode than the lower (1,1) mode, 

v) ^max slightly greater for a plate with K = 1 .9 than 
with X= 1.7, buckling in the same mode, under similar 
conditions. Figures 3.13(a) and (b) . 

In the present analysis, a viniform load has been applied 
whereas most of the available results in literature quote stress 
values for the case when a uniform displacement is applied on the 
edges. In the case of uniform load, seme small shear stresses do 
develop on the edges and hence the difference in the two results 
is not suJTprising. Actually, the uniform load at higher values 
(P » P _) / causes the loaded edge to bend inplane, with the 
central portion having moved longitudinally by a maximxam distance; 
this distance reduces towards the comer. Hence the stress values 
are likely to involve more error on this section (i.e. the loaded 
edge), than at other transverse sections. 



A 


Fig, 3.1 (a) Loaded plate with axes of symmetry 



Fig. 3.1(b) Discretization details of a quadrant . 
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Fig. 3.4 Load- deflection curves for an isotropic 
square plate, mode (1,1). 





edge of the plate centre of the plate 

(a) Symmetric configuration, mode (1, 1) 



(b) Anti -symmetric configuration, mode (2, 1) 

Fig. 3.5 Deflection profiles for an isotropic 
squre plate, FE condition. 
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edge of the plate centre of the plate 


(a) Symmetric configuration, mode(1,1) 



edge of the plate centre of the plate 

(b) Anti -symmetric configuration, mode (2,1) 

Fig. 3.6 Deflection profiles for an isotropic 
square plate, ME condition. 
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(b) A= 1.9 , mode (2,1) 

Fig. 3.7 Load -deflection curves for isotropic 





(b) FE condition 

Fig. 3.8 Load -shortening curves for 
an isotropic square plate. 






Fig. 3.9 Load -shortening curves for an 

isotropic square plate, rnoded, 1) 
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(b) A = 1.9, ME condition 
Fig. 3.10 Load-shortening curves for isotropic 
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(a) A =1.7^ FE condition 
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(b) A =1.9, FE condition 

Fig. 3.11 load-shortening curves for 
isotropic plates . 






Fig. 3.12 Effective width curves for isotropic 
square plates. 
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(b) A = 1.9, FE condition 
Fio 3 13 Effoctivo width curves for on isotropic 








(b) FE condition 

Fig. 3 14 Effective width curves for an 
isotropic square plate. 
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Fig. 3.15 Stress distribution in an isotropic 
square plate for load Nj( = 7.0, 
mode (1,1), ME condition. 
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Fig. 3.16 Stress distribution in on isotropic 
square plate for load Nx = 7.0, 
mode (1,1) FE condition. 
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Fia 317 Variation of maximum longitudinal 
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CHAPTER 4 


ORTHOTROPIC PLATES SUBJECTED TO LONGITUDINAL COMPRESSION 

4.1 INTRODUCTION 

In Chapter 3, the finite element method was applied 
to study the postbuckling behavior of isotropic plates. In 
this chapter, the scope of the finite element formulation has 
been enlarged to include the postbuckling behavior of speci- 
ally orthotropic plates. Such a plate or a laminate is obta- 
ined when orthotropic laminae are symmetrically arranged about 
the middle surface of the plate. 

The details of the analysis are similar to those for 
isotropic plates (Chapter 3, Section 3*1, Figure 3,1). The 
boundary conditions used are also the same# i.e. all the four 
edges are considered simply supported with unloaded edges 
either unrestrained (ME condition) or fully restrained (FE 
condition) against the normal (transverse) displacements. The 
relevant details of FEM have already been presented in Chapter 
2} the main difference lies in calculating the element stiff- 
ness matrices Wiich are quite involved for orthotropic plates. 
In order to study the effect of the degree of orthotropy 
E (= B /E ) on the postbuckling behavior# two materials are 

^ y ■ 

considered: glass-epoxy (E = 3.0) and boron-epoxy (E = 10.0). 

Results like critical loads# load-deflection curves, load- 
shor'lll^ning curves# changes in buckle patterns# effective width 



curves and stress distribution are reported. Wherever poss- 
ible, these results are compared with those of other investiga- 
tors . 


4,2 CURVES FOR BUCKLING LOADS 


For a simply supported rectangular orthotropic plate. 


subjected to uniaxial compression, the critical load N (per 

3CC XT 

unit width) is given by (Szilard, 1974) 


2 


N 


xcr 


= k 


U 7D D 
- — X..Z 


(4,1) 


where the compressive buckling coefficient k^ is defined as 
k^ = n^Td^Td^’ + 2H/V(D^Dy) (X%^) 

(4.2) 

in which (n = l), 

H = ■ + 2D^y (4.3) 

and = ^x’^y “ ^y^x (4.4) 

It is worth noting that the expression (4.2) is valid vAien the 

unloaded edges are either movable or are maintained straight 

(i.e. for ME and SE conditions only). The N (or k ) thus 

depends on the relative values of flexural rigidities D^, 

and shear rigidity The curves will, therefore, be 

different for materials having different degree of orthotropy 

E (= E /E or D /D ) , also called the rigidity ratio. For 
X y X y 

glass-epoxy (E =3) and boron-epoxy (E = 10), the curves are 



shown in Figures 4.1(a) and 4.2(a) for ME conditions. The 
absolute minimiam value of the critical load is obtained when 

1/4 

X = (4.5) 

y 

The transition from m to m+1 half waves occurs at 


\ = 'T ra(mtiy V(D^/Dy) 


(4.6) 


the successive values of X for glass-epoxy and boron-epoxy are 

1.86, 3.22, 4.56, and 2.51, 4.35, ..... respectively, for 

ME conditions . 

When the unloaded edges are rigidly hel|^ apart i.e. 
for FE condition, the previous expression for gets modified 
to (Harris, 1975) 

(mVx^) VTdTdT f 2H/v^(D D ) + (X%2) ^{dJd) 

■u ^ Jl V ^ 

1 + Vy(xVm^) 


(4,7) 


The curves obtained using Eq, (4.7), are shown in Figures 4.1(b) 
and 4.2(b). It can be seen that the transition from m = 1 to 
m = 2 and from m = 2 to m = 3 occurs at X « 2.2, 3.6 for glass- 
epoxy and atX= 2.6, 4.5 for boron-epoxy. The following addi- 
tional points may also be noted as a result of comparison of 
the curves for the two materials and for the two different 
inplane boundary conditions! 

(i) The statement that for a particular aspect ratio, 
is smaller for FE than for ME condition, which was 
valid for isotropic plates, is valid for orthotropic 
plates as well. 



(ii) For same inplane conditions on the unloaded edges 

and for the same aspect ratio, is smaller for boron- 
epoxy than for glass-epoxy; thus k^ decreases with the 
increase in the degree of orthotropy. 

(iii) The absolute minimrim value of k^, for a fixed m, 

occurs at a larger value of \ for materials having 
greater value of E. For example for m = 1, X is 1.316 
for glass-epoxy and is 1.778 for boron-epoxy under ME 
conditions. The corresponding values of k^ are 3.430 
and 2.605, Similar trend is also seen for FE condition, 

(iv) The transition from one curve to another occurs at 

different points (X value) for different inplane boun- 
dary conditions; the difference is more for weak 
orthotropic material (glass-epoxy) and less for strong 
orthotropic material (boron- epoxy) , The difference is 
maximum for isotropic plates, (Figures 3,2 and 3.3). 

(v) Higher the ratio E, smaller is the critical load for 
plates of same X' (= m/X) , under similar inplane 
boundary conditions. For E > 1, is less than 

that for E = 1 (i.e. isotropic plate). 

(vi) A plate with, X = 2,0, and E = 3.0 (glass-epoxy) 

buckles in two half waves for ME condition but in a 
single half wave for FE condition. ,The same plate 
with E = 10 (boron-epoxy) bucxles in a single half 
wave only under either of the inplane boundary condi- 
tions on unloaded edges. 



4.3 ILLUSTRATIVE EXAMPLES 

The two materials considered are glass-epoxy and boron- 
epoxy whose properties are: (Agarwal and Broutman, 1980) 

Glass-Epoxy 

E^ = 42 GPa, = 14 GPa, G^^ = 7 GPa, v = 0.25, v = 0.0833 

Boron- Epoxy 

E = 210 GPa, E = 21 GPa, G = 7 GPa, v = 3.0, V = 0.03 

#A- j ^ j 

For each material, the following illustratiye examples 
are taken: 

(i) Square plate (X=1.0), l.Omxl.Om and 10 mm thick — 
with unloaded edges free to move (ME condition), buck- 
ling in (1,1) and (2,1) modes. 

(ii) Rectangular plate (k = 2.0), 2.0 m x 1.0 m and 10 mm 
thick 

(a) unloaded edges rigidly held apart (FE condition) , 
buckling in vi,l) and (2,1) modes 

(b) unloaded edges free to move (ME condition) , buck- 
ling in (1,1) and (2,1) modes (for boron-epoxy 
only) . 

(iii) Rectangular plate (X = 3.0), 3.0 m x 1.0 m and 10 mm 
thick — with FE condition, buckling in (2,1) and (3,1) 


modes. 



4 . 4 NUMERICAL RESULTS 


4.4.1 Critical Loads 

For each illustrative example considered# the critical 
load is obtained from the load-deflection/shortening curves as 
explained in Section 3.3.1. The values of the buckling coeff- 
icient k^ are tabulated in Table 4.1 therein the exact values# 
as obtained from Eqs. (4.2) and (4.7), are also presented. 

The two sets of value are in very close agreement, 

4.4.2 Load-Def lection/Shortening Curves 

The load-deflection curves# between the load N 

2 2 — — 

(= N b /it Vd D ) and the out of plane deflection w (>= w 
X X o o 

Wq being the maximxim deflection) are drawn in Figures 4.3-4, 5, 
The curves in Figures 4.3(a) and 4,3 (b) are for square plates 
of glass-epoxy and boron-epoxy respectively# buckling in the 
primary mode (1#1) under ME conditions. The deflection pro- 
files for the symmetric and the anti- symmetric configurations# 
corresponding to modes (1#1) and (2,1)# are given separately 
in Figures 4.6 (a#b). In order to make a direct comparison of 
the results# the load-deflection curves are presented in 
Figure 4,7 between the load ratio N (= N /n or * N /k ) and 
w^. The results for the corresponding isotropic case are also 
reproduced. When compared with the results of Prabhakara arid 
Chia (1973) , some discrepancy is observed between the two 
results (Figure 4, 7 a!)» the deflections in the present analysis 
are on the lower side by about 15-20% # under ME conditions. 

The reasons for this difference are manifold: first 



4^1 Comparison of k values for orthotropic plates 
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Prabhakara and Chia (1973) have not specifically mentioned the 
type of inplane boundary condition used. Secondly, they have 
applied the condition of zero shear on the edges. Thirdly, 
the material properties used by them slightly differ from 
those used in this work. 

The load- shortening curves are drawn between the non- 
dimensional load and the non-dimensional shortening A 
(= AE b t/a % D ) and are shown in Figures 4.8-4.12. 

Figure 4.9 shows such curves for glass-epoxy plates with X = 

2.0 and 3.0 under FE conditions. These curves match with those 
of Chandra and Raju (1973). The small difference is because of 
different inplane conditions. Figure 4.12 gives a comparison 
of load- shortening curves for three materials, drawn on the 
same plot, for a fixed value of X' = 1 (X* = m/X). Again the 
trend of these curves is very much similar to that reported by 
Chandra and Raju (1973). 

4.4.3 Change in Buckle Pattern 

The possibility of a change in the buckling pattern 
is indicated by the point of intersection of the corresponding 
load- shortening curves. Whereas a clear point of intersection 
is obtained for a glass- epoxy square plate xmder ME condition 
(Figure 4.8 a ), there is no such phenomenon in the case of 
boron-epoxy plate (Figure 4.8 b). Figure 4,9 shows results 
of glass-epoxy rectangular plates of aspect ratios 2.0 and 3.0 
under FE conditions. In both cases, the changeover points 
[for X = 2.0 from mode (1,1) to (2,1) and for X = 3,0 from 
(2,1) to (3,1)] are distinct. Similar curves for boron- epoxy 

eo-ttwmenXS JS4*, 



plates are shown in Figures 4,10 and 4.11; the later figure 
shows divergence of the two lines corresponding to Cl,l) and 
(2,1) modes. 


4.4,4 Effective Width Curves 

The effective width formula for isotropic plates, 
expressed by Eq, (3,9), gets modified to (Chandra and Raju, 1973) 


b 

e 


N a 

X 

At E^ 


(4.8) 


for orthotropic plates. Curves are drawn to represent varia- 
tion of non-dimensional effective width b with the non-dimen- 

e 

sional load as well as with the load ratio N, Figure 4,13 
is for glass-epoxy plates of X = 2.0 and 3.0 under FE conditions 
and Figure 4,14 for a boron-epoxy plate of X = 2.0 with ME and 
FE conditions. The variation of effective width with the para- 
meter X' (= ^/pt) is illustrated in Figures 4.15 and 4.16. For 
a fixed value of X' (« 1.0), Figure 4,15 shows the variation 
under different inplane boxindary conditions idiereas for a given 
inplane boundary condition (FE) , the variation of for 

different values of X' is shown in Figure 4.16, Results of 
Figure 4, 15(b), (X' = 1.0, FE condition) are quite close to 

those of Chandra and Raju (l973) . 

A 

4.4,5 Stresses 

For the problems analysed, stresses have also been 
calculated at some load levels. The stress distribution 
follows a similar pattern as in Figures 3.15 and 3.16 for 



isotropic plates / \inder the two inplane boundary conditions. 

Figures 4,17-4.19 show the variation of non-dimensional maximum 

— 2 2 

longitudinal stress d (= d b t/% DD) with the non- 

irioL^ HtcLX. X y 

dimensional load N . The variation in postbuckling range is 
also linear, as in an isotropic plate. 

4.5 DISCUSSIONS 

4.5.1 Load-Deflection/Shortening Curves 

From the load-deflection cuirves in Figures 4. 3-4. 7, 
following conclusions are drawn! 

(i) For the same load ratio N, the deflection w^ is 

higher for higher values of the rigidity ratio E of 
the plate (Figure 4,7). It is obvious from these 
curves that, other parameters remaining same, the 
deflection of an orthotropic plate is greater than 
that of an isotropic plate. 

(ii) The deflection curve for the glass-epoxy plate (E = 

3.0) is only marginally different from that of an iso- 
tropic plate for ME conditions and the two curves 
almost coincide for FE conditions (Figures 4.7 c and 
4.7 d . This fact is also obvious from Figure 4, 12(b), 
where N - A lines for the isotropic and the glass- 
epoxy plates for the same are parallel. However 
the curves of boron-epoxy are quite distinguished from 
other materials, under both inplane conditions. 
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(iii) For the same X' (= m/X) , the deflections are larger 
for ME conditions than those for FE conditions 
(Figures 4.7a, d) . This is also clear from Figure 4.4, 
in which N - w curves are drawn for boron-epoxy with 
X = 2.0, under ME and FE conditions. 

(iv) Whereas for a plate of a fixed aspect ratio, under 
ME conditions, the tendency for the maximum lateral 
deflection is to decrease from (1,1) to (2,1) mode 
(Figures 4.7a, b) , the curves of Figures 4.7c, d (for 
X = 2.0, under FE conditions) show just the opposite 
trend, whereby the maximum deflection in (2,1) mode is 
greater than that in (1,1) mode. 

Further^ from the load- shortening curves in Figures 
4.8, 4.10 and 4.12, it is observed that: 

(v) For a plate of a particular material, size and mode 
the reduction in the axial stiffness, after primary- 
buckling, is larger for ME condition than for FE 
condition (Figures 4.10a, b and 4.12a, b) . 

(vi) The reduction in -the axial stiffness increases as the 
value of E increases (Figures 4.8a, b and 4.12a, b) . 

4.5.2 Change in Buckle Pattern 

Following points regarding change in buckle pattern, 
are observed with reference to the various load- shortening 
curves shown in Figures 4.8-4.11: 

(i) In the case of boron-epoxy with X = 2.0 and under FE 

condition , the load- shortening lines in (1,1) and (2,1) 
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modes, (Figure 4.10b) intersect at the load N = 4.43, 
a value slightly higher than the critical load corres- 
ponding to the secondary mode (2,1), viz. 4.00. On the 
other hand under ME condition, (Figure 4.10a), the 
point of intersection is at 6.40, which is much above 
the critical load of (2,1) mode, viz. 4.20. Similarly 
for glass-epoxy plates having X = 2. 0 and 3.0 (Figures 
4.9a, b) the possible changeover appears to be at a load 
slightly greater than the critical load corresponding 
to the secondary mode. Thus under FE condition, the 
change in buckle pattern is rather quick (as observed 
for isotropic plates also, in Chapter 3). 

(ii) Figure 4.11 is the graph of N against (A/A^^.) for a 

boron-epoxy plate with X = 3.0. It is clear that (1,1) 
and (2,1) lines do not meet, implying that the change- 
over from a higher (2,1) to a lower (1,1) mode is not 
possible. 

(iii) The approximate formula to predict the point of inter- 
section given for isotropic plates, Eq, (3,15), seems 
to hold good for orthotropic plates, also. However a 
more rational formula, which takes into account the 
parameter m as well, is proposed and is valid for ME 
condition: 
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where and are, respectively, the critical 

loads corresponding to the secondary and the primary 
modes, N . is the load point at which the two N - A 
lines corresponding to the primary mode (= m) and the 
secondary mode C= intersect. The value of m, for 

example, will be 1, for changeover from Cl,l) to (2,1) 
mode and will be 2 for changeover from (2,1) to (3,1) 
mode . 

An interesting observation is made for - A lines 
of boron-epoxy square plate under ME condition (Figure 4.8b). 
The two lines, as per the prediction rule described in Section 
3.4.2, should have met whereas they appear to be almost para- 
llel. Actually the two lines do meet but at a very high load 
which is not of any practical significance# 

4.5.3 Effective-Width Curves 

Empirical formulae to obtain effective width for ortho- 
tropic plates are not many. As in the case of isotropic plate, 
the formulae are based on the assximption that the mode at the 
instant of buckling remains unchanged. However in the present 
investigation, effective width curves, presented in Figures 
4.13 and 4.14 incorporate the possibility of the change in 
buckling pattern. 

Following important features can be noted from these 
curves and those of Figures 4,15 and 4.16: 

(i) The changeover from a lower to the next higher mode is 
quicker vinder FE condition than mder ME condition 
(Figures 4,14a, b). 



(ii) For a particular load ratio N, the effective width 

reduces with the increase in the value of X' (Figures 
4.16a, b and 4.15 b) . 

(iii) For the same X' and under same inplane conditions 

the effective width b decreases with the increase in 

e 

the rigidity ratio E (Figures 4.15 and 4.16). 

(iv) The reduction in effective width is more in the 

initial post-buckling range and is very small for the 
load ratio N > 1.75 (Figure 4.15). A similar trend 
was also observed for isotropic plates (Figures 3.12- 
3.14). 

4.5.4 Stresses 

Following are some of the general remarks made on 

various stress results obtained in this analysis. (These also 

apply to the isotropic plates) J 

(i) The maximum value of membrane (longitudinal) stress, 

a , occurs- at the crest of the buckle on the 
max 

unloaded edges, under ME condition and at the corners 
of the loaded edges under FE condition. 

(ii) The difference between the values of at the crest 

of the buckle and at the edge comer, for a particular 
load level is very small under FE condition but, is 
quite large under ME condition. 

(iii) Under the FE condition, when the plate buckles in (l,l) 
mode, the transverse (membrane) stress develops on 



90 


the THiloaded edges J this stress is con^ressive near 
the sides and tensile in the central part of the edge. 
Some additional remarks are made with reference to 
Figures 4.17-4.19: 

(iv) The variation of d with the load N is linear in 

max X 

the postbuckling range as has been the case with 

isotropic plates also (Figures 4.17-4.19). 

(v) The value of o relative to the value at the instant 

max 

of buckling 5^^ (. t/k^ '^5;^). reduces with 

the increase in degree of orthotropy E. This is 

clearly seen in Figures 4.17a# b, vhere the N “ a 

line for boron-epoxy (E = 10) is having greater slope, 

than that for glass-epoxy (E = 3) for the same mode. 

(vi) When graphs of Figures 4.18 a and b are compared 

for the same mode (2,1), it is found that, for a given 

load ratio, the maximum stress ratio Cf /d , (or 

d /d ) reduces with the increase in the aspect ratio, 
max: cr 

(vii) From Figure 4.19, it can be interpreted that for a 

given plate, buckling in a particular mode, the value 

of d , is greater for ME condition than for FE cond- 
max 

ition, at a specified load ratio. 




(a) ME condition 



Fig. 4.1 Compressive buckling coefficients 
for glass -epoxy plates! 
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Fig. 4.2 Compressive buckling coefficients 
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(a) Glass-Epoxy, FE condition, mode(2, 1) 



(b) Boron-Epoxy , FE condition, mode (2, 1) 
Fig. 4.5 Load - def lecti on curves. 
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Fig. 4.6 Deflection profile for a square plate, 
ME condition. 
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B- Boron-Epoxy (E = 10) 


Fig. 4.7 Load -deflection curves between load ratio 
and maximum deflection. 
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Fig. 4.9 Load -shortening curves for glass 
epoxy plates. 
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Fig. 4.10 Load “Shortening curves for boron 
epoxy plate. 
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Fig. 4.13 Effective with curves for glass 
epoxy plates. 
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Fig. 4.14 Effective width curves for a 
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Fig. 4.18 Variation of maximum longitudinal 
stress with the load for glass- 
epoxy plates . 
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Fig. 4.19 Variation of maximum longitudinal stress 
with the load for boron-epoxy plate. 







CHAPTER 5 


PLATES SUBJECTED TO SHEAR LOADING 

5.1 INTRODUCTION 

For buckling to take place# it is not necessary that 

the plate be loaded in compression only, A system of forces 

that can cause compressive stresses anywhere in the plate is 

responsible for its instability. The pure shear# for example# 

o 

produces compression on planes inclined at 45 and when these 
forces become increasingly large# buckling takes place. Since 
the compressive stresses are developed on planes that are 
inclined# the buckle pattern tends to become skew. In case of 
long plates# it is reported (Gerard# 1962) that the half wave 
length equals to 1.25 times the width. 

From the review of literature# presented in Chapter 1, 
it is seen that not much work has been done on the postbuck- 
ling analysis of shear loaded plates. Moreover# the Finite 
Element Method has hardly been attempted. In this . chapter# 
simply supported isotropic and orthotropic rectangular plates# 
subjected to pure shear load# are analysed with the view to 
study their postbuckling behavior. 

5.2 ANALYSIS DETAILS 

In Chapters 3 and 4# only a quadrant of the plate yra.s 
considered because of the syrametiry with respect to loading as 
well as geometry. In case of a plate subjected to pure shear, 
symmetry with respect to loading is violated. Hence, the full 
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plate needs to be analysed, as shown in Figures 5.1. It is 
discretized in 16 equal 4-noded rectangular finite elements. A 
uniform load (in x~direction) along the edge y = b (in the form 
of concentrated loads at the edge nodes) is applied, while the 
edge y = 0 is restrained to move in its own direction, thus 
causing a pure shear. The other pair of edges viz. x = 0, a are 
not allowed to- translate in the y-direction. These inplane 
conditions correspond to the NASTRAN model as reported by 
Agarwal (l9Sl). The plate is assumed simply supported on all 
the four edges and hence the appropriate boundary conditions of 
zero deflection and zero moment, namely the Eq. (1.4) are applied.; 

Buckling Pattern/Mode 

Unlike a longitudinally compressed plate, a plate subjects 
to shear loading, does not buckle in a simple (uncoupled or couple 
mode with n equal to unity. The buckling pattern is actually a 
mixtvire of various uncoupled modes with different combination of n 
and n values. In the case of an isotropic plate, it is well know: 
that the resulting buckled pattern is such that the sum (m+n), fo| 
each constituent (uncoupled) mode is either even or odd. When (ml 
'is even, it requires the slopes to be zero ( 3w/ 9x = 3w/3y = O) at; 
the centre of the plate, while for odd values of (m+n), the condi 
of zero deflection Cw = O) is needed. The two patterns (modes) I 
may, therefore, be called as symmetric and anti -symmetric, respec 
ively. Hence, for a plate of given aspect ratio, either of thasa 
modes will yield the minimum critical load. | 
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To initiate a synunetric mode, a small lateral load is 
applied at the centre, '0' of the plate _( Figure 5.1a) and for 
an anti-symmetric mode, two loads (in opposite directions) at 
nodes M and N are applied as has been the case with plates under 
edge compression. But, the condition of zero slope for the 
symmetric mode and the zero deflection condition for the anti- 
symmetric mode are imposed only at the centre of the plate 0 and 
not on the whole line PR, so that buchling pattern is free to 
take its own shape. 

5.3 CURVES FOR SHEAR BUCKLING COEFFICIENTS 


For a simply supported, rectangular, isotropic plate, 
under the edge shear, the critical load (per unit length) 

is given by (Kumar, 1985) 


where k^ is 
shear stress 


N = k 

xycr s 


rr^D 

72“ 


(5.1) 


the shear buckling coefficient. In terms of critical) 


T the above formula may be re-written as : 
cr 


Ji^E 

T = k„ — — 

12(1 - v) 




(5.2) 


In these equations b is always the shorter dimension of the j 
plate, as all the edges are loaded in shear. Hencd, in this casi 
the aspect ratio X(= a/b) < 1 does not have any meaning; X will) 
always be > 1*0. For an orthotropic plate, the corresponding | 
expression for the critical load is 



( 5 » 3 ) 


N 

xycr 


71 n/'D D 

k ^2L Z. 


In Eq. (5.1) or (5.3), the expression for k_ is not a simple 
function of m, n, and x » unlike the case for plates vinder edge 
compression; in fact has to be evaluated numerically for each 
case. 


Based on Eq. (5.1), curves for the shear buckling coefficien 
kg, for isotropic plates, are shown in Figure 5.2. 

Following important points are worth noting with respect 
to these curves : 

i) Curves start from X = 1 *0 and the value of k is highest 
at this aspect ratio. The minimum value of k_ is 5.35 

O 

for X -» (i.e. a long plate) . 

ii) Each curve is the combination of curves corresponding to 
m (= 1,2,3, •»•) and n (= 1,2,3, •••) number of half waves 
in X- and y-direction, respectively, 

iii) The s^munetric pattern (of one large buckle) and the anti- 
symmetric pattern (of two large buckles) intersect at 
X = 2.05 and thus frcxn X = 1 .0 to 2.05, the symmetric | 

pattern gives the minimum critical load. Similarly, frcm | 
X= 2.05 to 3.50 (where the next two curves of anti-symmetrj 
and symmetric patterns meet), the anti-symmetric pattern 
governs the buckling, and so on. I 

The critical load and even the postbuckling path is j 

independent of the direction of the shear (Figures 5.1b,c) for 
isotropic as well as for specially orthotropic plates. However, i 
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for anisotropic plates (not discussed here) there are two 

different critical loads and the postbuckling paths for the 

two alternative directions of the shear (Zhang and Matthews, 1934) • 

5.4 ILLUSTRATIVE EXAMPLES 

Following Illustrative examples are taken in order to 
study the postbuckling behaviour of isotropic and orthotropic 
plates under shear : 

I : scruare plate (X= 1 «0) / 1.0 m x 1.0 m and 10 mm thick, 

E = 200 GPa , V = .316, buckling in symmetric and anti- 
symmetric modes. 

II ; rectangular plate (X= 1*5) / 1.5 m x 1.0 m and 10 mm 

thick, E = 200 GPa , V = .316 , buckling in symmetric 
and anti -symmetric modes- 

III i rectangular plate (X= 2-5 )/ 2-5 m x 1-0 m and 10 mm 

thick, E = 200 GPa , V = -316, buckling in symmetric 
and anti -symmetric modes. 

IV : Glass-Epoxy square plates (X = l.O) / 1.0m x 1 .0 m and 

10 mm thick, E = 42 GPa, E = 14 GPa, G = 7 GPa, 

X y xy 

V = .2 5, V = .0333, buckling in symmetric and anti- | 

^ y I 

symmetric modes. 

V : Boron-Epoxy, sguare plate (X = 1 .0) / 1.0 m x 1 -0 m and 

10 mm thick, E = 210 GPa, E = 21 GPa, G = 7 GPa, 

X y xy I 

^x “ *20, Vy = .03, buckling in symmetric and anti- 
symmetric modes - 
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5*5 NUMERICAL RESULTS 
5.5«1 Critical Loads 

As before, the critical loads are obtained from the load- 
deflection/shear displacement curves (ref. Section 3.4.1), for 
each example. The values of shear buckling coefficients k^ are 
tabulated in Table 5.1. These values for an isotropic plate, 
are compared with those obtained by Stein (1947), both for 
symmetric as well as for anti -symmetric modes. The agreement, 
between the two set of values is found to be very good. 


Table 5.1 Comparison of shear buckling coefficients k^ 
Isotropic plates 


X 

Mode 

, Value of k^ 

3 

Exact 

( Stein, 1947) 

Present 

work 

1 .0 

S 

9.35 

9 .43 

1 .0 

A 

11 .63 

11.76 

1.5 

S 

7.07 

7 .11 

1.5 

A 

7.97 

S .04 

2 .5 

S 

6.29 

6.34 

2 .5 

A 

6 .06 

6.12 


S : Symmetric A : Antisymmetric 

O rthotropic plates 

\ = 1.0, Glass-Epoxy, Symmetric Mode : k_ = 8.65 
X =1.0, Eoron-Epoxy, Symmetric Mode : k^ = 


9.05 
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5.5*2 Load-Deflection/Shear DisplacCTient Curves 

The applied shear load N along the edges is non- 

xy 

dirnensionalized as = N b^/n^D for isotropic plates and 

b /n Dy for orthotropie plates. The load is also 

expressed as the load ratio N (= N /N = N /k ) , where 

xy xycr xy s ' 

^xycr critical shear load. Load-deflection curves between 

N and (= w^/t; w^ being the maximum lateral deflection, as 

usual) are shown in Figure 5.3, for isotropic plates, having 

X = 1.0, 1.5 and 2.5. Both, symmetric as well as anti -symmetric 

modes are considered. Figure 5.4(a) shows such curves for glass- 

epoxy and boron-epoxy square plates. To have an overall picture 

of the shape of the buckle patterns, the deflection contours are 

drawn on the plan of the plate, (Figures 5.6 - 5.S) . 

Corresponding to the term 'shortening' ( A), in a 

longitudinally compressed plate, a term 'shear displacement' 

( Agj^) , has been used in the present case of a shear loaded plate. 

This may be defined as the relative displacement of the longitu- ; 

dinal edges, parallel to itself, and is normalized as A^^ 

(= A ^ G bt/r D) for an isotropic plate and as i 

s n 

(A , G b t)/(7r ''D D ) for an orthotropic plate. Plots of 

sh xy 5C y 

load versus shear displacement ^sh are presented in Figure j 
5.5 for isotropic and in Figure 5.4(c) for orthotropic plates 
(glass-epoxy and boron-epoxy square plates buckling in symmetric; 
modes). ^xy ” ^sh curve for a plate with X = 2.5 (Figure 5. i 

is close enough to that of Stein (1985), who has analysed a longi 
plate only with edges held_ straight . 
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5.5.3 Stresses 


For each illustrative example, stresses have also been 
calculated. The shear stress values have been marked on the 
edges of the plates in Figures 5.6 - 5.3. Figure 5.9 depicts 
the variation of maximxim shear stress with the load, for plates 
with X = 1.0, 1.5 and 2.5 and buckling in symmetric or anti- 
symmetric modes. Similar ciirves for glass-epoxy and boron-epoxy 
square plates are shown in Figure 5.4(b). 


The direct membrane stresses a and a are zero, in a 

X y 

shear loaded plate, until buckling starts. In the postbuckling 
range, these stresses, which could be tensile or compressive, 
also get developed. The maximum values of longitudinal tensile 
stress cr^ max transverse tensile stress are normalized: 

with respect to the critical load (of uniaxial compression) 

“x max = ‘’x max '"■'“xcr °y max = ‘’y max ^/“xcr’ respective!; 
These are plotted against the load ratio N in Figure 5 #10, for 

isotropic, glass-epoxy and boron-epoxy square plates, buckling ip 
symmetric mode. The magnitude of the longitudinal tensile stress 
a , produced in an isotropic plate with X = 2.5, at the load 
N = 1.5 is canpared with the long plate results of Stein (1935). | 
This value in the present case is 0.39 v^ile for the long 

plate it is approximately 0.30 The difference in the two ' | 

values is, perhaps, due to the fact that Stein maintained the | 
loaded edges of the long plate to be straight. i 
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5.6 DISCUSSIONS 
5 .6 .1 Deflections 

Some conclusions drawn from the various deflection curves/ 
for the shear loaded plates are : 

(i) The deflection increases with the aspect ratio 

of the plate for the same load ratio N and for the same 
mode (Figure 5.3). 

(ii) For a particular value of N the deflection w^ is less 
for a shear loaded plate than for a plate loaded in 
compression (Figures 5.3a . 4.7a). 

(iii) In regions near the corners of the compression diagonal i 
of a buckled plate the lateral deflections are of 
opposite nature to that of the rest of the plate. 

(iv) For the loads considered (upto about twice the critical 
load value)/ it is observed that the load-shear 
displacement lines, for all aspect ratios, are almost 
parallel (Figure 5.5). 

(v) There is only a slight deviation of the postbuckled i 
line from the original (prebuckled) line. This shows 
that for a shear loaded plate, the reduction in the 
shear stiffness after buckling is very small (about 119^) 

(vi) As the N - A. lines (for loads upto twice the critical 
xy ai I 

corresponding to the two modes (symmetric and anti- I 
symmetric) do not meet for the plate of any aspect ratic 
(Figure 5.5); the question of change of wave pattern fri 
a symmetric to an anti -symmetric mode (or vice versa) d< 
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not arise ’Stein (19S5) also has stated that the mode 
at the instant of buckling remains same, in the initial 
postbuckling regime. 

(vii) For orthotropic plates (glass-epoxy and boron-epoxy) 
the deflections are greater than those of an 
isotropic plate for the same load ratio N and the mode; 

increases with degree of orthotropy S, (Figures 5.3a, 
5 a ) • 

(viii) The reduction in shear stiffness decreases with the 

increase in value of E/ thus the shear displacement for 
a boron-epoxy plate is less than that for a glass-epoxy 
and for an isotropic plate, for the same load N 

xy 

(Figures 5.4c, 5.5). It may be noted here that just the 
opposite observation holds for orthotropic plates 
loaded in compression. 

Stresses 

Following points are observed frctn the results on stresses r 

The variation in shear stress "C is symmetrical with 

xy 

respect to the two diagonals of the plate. 

The maximum value of the shear stress occurs at the 

two corners (diagonally opposite) , for isotropic as well 

orthotropic plates, buckling in the mode corresponding to- j 

the minimum critical shear load, (e.g. for X= 1»0, 1.5 in; 

symmetric mode and for X = 2.5 in anti-symmetric mode). 

T is minimum at the other two comers, 
xy I 


5 .6 .2 

(i) 

(ii) 
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(iii) 


(iv) 


(v) 


( Vi) 


(vii) 


For plates of low aspect ratio (X = 1.0, 1.5), buckling 
in an anti -symmetric mode, T occurs at some distance 

itla X 

(about one quarter) from the above comers (Figures 5.6, 
5.7) . Thus its position is at nodes 2,6,20 and 24 
( Figure 5 .la) . 

The magnitude of the maximum longitudinal tensile stress 

O’ and transverse tensile stress O increases 

X max ■ y iinax 

linearly with the load (Figure 5.10). For a particular 

value of the load ratio N, O increases and a , 

x max y max 

decreases with the increase in the value of E 
(Figure 5.10a,b). 

For the same value of N, the stresses 0^ max niax 

reduce as the aspect ratio increases. 

For a fixed value of the shear load N , the maximum 

xy 

shear stress X is areater in symmetric mode than ' 
max 

in anti— symmetric mode, for plates with X = 1.0, 1*5 and 
it is greater in anti -symmetric mode than in symmetric 
mode, for the plate with X = 2 .5 (Figures 5.6 - 5.S) . 
This shows that the magnitude of the absolute maximum 

shear stress occurs when a plate buckles in its funda- j 

! 

m.ental mode (mode in which the critical load is absolute 
minimum. | 

From Figure 5.8, it can be interpreted that for a given! 
load ratio ii, the magnitude of x decreases with the; 

niaX 

increase in the aspect ratio. For example, the values ; 

X for X= 1.0, 1.5 and 2.5 are approximately 27.0, 20 
max i 

and 14.4 respectively. 



121 


C--*Pj 


V =0 on BC 


P, =Nxy(a/4) 
^2 = 2 


u = 0 on AB 


B X 


Fig. 5.1 (a) Rectangular plate under 
shear loading. 



Fig. 5.1(b) Positive Fig.5,1(c) Negative 
shear . shear. 




X (=a/b) 

Shear buckling coefficients for simply supported 
isotropic plates. 
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A = 1.5, (c) Symmetric 
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(b) Anti-symmetric 
mode 
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(d) Anti-symmetric 
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Fig. 5.3 Load-deflection curves for isotropic 
plates loaded in shear. 








curves shear stress with load 



Fig. 5.4 Various curves for orthotropic 
square plates , loaded in shear, 
buckling in symmetric mode. 
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^sh 

Fig. 5.5 Load -shear displacement curves, for 
isotropic plates, loaded in shear. 
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Fig. 5.6 Deflection contours and shear stress distri- 
bution for isotropic square plate loaded in 
shear. 








Fig. 5.7 Deflection 
stress di; 


[tllM 
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I contours and shear 
stribution for isotropic 
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Fig. 5.9 Variation of maximum shear stress 
with the load . 
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(a) Longitudinal stress 


max 


(b) Transverse stress 

Fig. 5.10 Variation of maximum tensile stress 
with the shear load for a square 
plate buckling in symmetric mode . 






CHAPTER 6 


SUMMARY AND CONCLUSIONS 


6 .1 SUMMARY 

The main aim of the present work has been to investigate 
the postbuckling behavior of isotropic and orthotropic thin 
rectangular plates subjected to either longitudinal compression 
or edge shear* Plates are assximed to be simply supported on 
all edges* Two inplane boundary conditions on unloaded edges 
have been considered : (i) edges free to move (ii) edges fully 

restrained. Throughout the investigation, the finite element 
method has been employed/ a 16 element mesh is used to discretiz* 
the quadrant of the plate (for compression loading) and the full! 
plate (for shear loading)* Each element is a 4-noded rectangle | 
with 6 degrees of freedom at each node. The governing different 
equations have been derived by minimizing the total potential 
energy of the element and then assembled to obtain a set of 
simultaneous nonlinear algebraic equations which have been solve| 
by an incremental-cum-iterative method. 

■ i 

Main contributions of the present study are : 

(1) It has been shown that the finite element method is quite! 

i 

versatile in handling complex non-linear problems sucdi as; 
one dealt in this thesis. 

( 2 ) A method has been suggested to predict a possible change ! 
in the buckling mode for a given size of the plate and | 
for specified boundary conditions. 
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C3) Empirical formulae have been deviced to predict the 
magnitude of the load at which a possible driangeover 
in the buckling mode takes place. Equations (3»15) and 
(4*9) are applicable for longitudinally compressed plates 
and are valid for isotropic as well as orthotropic plates 
and for any inplane boundary condition# 

( 4 ) A useful set of effective width curves has been presented 
these curves are valid even after the change in the 
buckled configuration has taken place# 

(5) The FE inplane boundary condition has been explored in 
detail for the postbuckling analysis; not much work is 
available in literature for this condition# 

(6) By investigating plates of different aspect ratio# under 
shear loading# it has been shown that the primary mode 
at the instant of buckling remains unchanged in the 
postbuckling range# The buckling mode is a coupled one 
consisting of either symmetric uncoupled modes or anti- 
symmetric uncoupled modes# 

6.2 CONCLUSIONS 

Following useful conclusions can be drawn for plates 

subjected to longitudinal compression : 

(i) The values of the critical load obtained using finite 
element method are very close to the / exict o nesj ^ 

Cii) The postbuckling behavior is greatly influenced by 
inplane bovindary conditions on unloaded edges# 
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(iii) 


( iv) 


(v) 


(vl) 


(vii) 


A plate of given aspect ratio may buckle in different 
modes under different inplane conditions. For example, 
an isotropic plate with X = 2 .0 buckles in (2,1 ) mode 
under ME condition but in (l,l) mode under FE condition 
A plate initiated to buckle in a particular symmetric 
mode, converges, for a given load, to a mode which is 
a symmetric but not necessarily to the same one. On 
the other hand, a plate initiated to buckle in a 
particular anti -symmetric mode would finally converge 
to the same anti-symmetric mode. 

For a plate of given aspect ratio X and a specified 
mode (m,n), i.e. for a fixed value of X' (= m/X), the 
maximum lateral deflection is greater under ME condition 
than that under FE condition. 

For the same load ratio, the lateral deflection of an 
orthotropic plate is greater than that of an isotropic 
plate and it increases with the value of rigidity ratio 
E of the material. 

The axial stiffness of a plate after buckling reduces 
to about to lCf)i of the prebuckled value, the actual ; 
value being dependent on plate boundary conditions, plate 
dimensions and the mode of buckling. A f\irther reductiorj 
takes place after the secondary buckling to about 4096 ; 

to 5096 of the original (unbuckled) value. The reduction | 
in stiffness is greater under ME condition than under 
FE condition . 
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(viii) Ihe errective vd.dth reduces with an increase in the 
axial load; the reduction is large enough in the 
beginning but is very small for loads beyond 1.75 P 

cr 

(ix) The efrective vddth of a plate is greater under FE 

condition than under ME condition/ heeping other factors 
same. 

(x) For a given X' and for prescribed boxmdary conditions, 
the effective vd.dth b^ decreases with an increase in 
the degree of orthotropy E. 

(xi) In tlie postbuckling regime/ a plate undergoes differen- 
tial shortenings causing non-uniform strains and 
stresses. The longitudinal stress under ME condition 
is maximiom at the ere so of the buckle whereas, for FE 
condition/ it is raaximum at the comers (junctions of 
loaded and unloaded edges) . 

(xii) For the same load ratio, , decreases with the incre- 

ase in the value of E. For example, the ratio G" /cr~ 

'^max'^cr 

for boron-epoxy is lower than that for glass-epoxy. 

Additional comments on ' chanqe in buckle pattern* of the plate 
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was then obtained as point of intersection of the corresponding 
load— shortening lines. However, Supple (1966) has found that 
the change in buckle pattern is in the form of a coupled mode 
and it branches from the secondary mode (i.e. the one with 
higher critical load) . But, for plates under FE condition. 
Supple (1970) shows that rhe change from primary to secondary 
mode does taJce place. 

Since the present work does not examine coupled modes, 
the author feels that the results of effective width curves of 
Chapters 3 and 4 (involving change in buckle pattern) are 
approximate only. The curves however represenr the qualitative 
picture as no how these curves get effected if the change in 
buckling pattern is to be accounted for. 

Following useful conclusions can be drawn for plates 
sub j ected to edge shear : 

(i) The mode at the instant of buckling is a mixed one 
consisting of eitlier uncoupled symmetric modes or 
uncoupled anti- symmetric modes. 

(ii) The load-shear displacement lines for a plate of a 

given X are parallel; hence a change in buckled pattern 
is not possible, atleast in the postbuckling range. 
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investigated (p < 2 . This has been found true 

for all the three values of X considered. 

(ill; The reduction in the shear stiffness after primary 
buckling is rather small. 

Civ) As expected, the deflection profiles are skewed 
making seme angle with the edges. 

(v) Longitudinal and transverse tensile membrane stresses 
are produced which tend to stiffen the plate in its 
postbuckling range. 

( vi ) Distribution of shear stress is -symmetrical with 
respect to the two diagonals of the plate, with 
maximum value occuring at the two ends of the elongated 
(tension) diagonal. 

6.3 SUGGESTIONS FOR FUTURE WORK 

There is a tremendous scope of further research 
the subject of postbuckling behavior of plates. Some sug 
in this regard are ; 

# The effect of different flexural boundary cond 
should be studied, both for isotropic and orthotropic pla 

* The postbuckling behavior of laminated composi 
plates should be studied. 

if The effect of combined loading should be inves 
so as to make the problem more realistic. 
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* The effect of transverse shear in the plate theory 
can be included to make the large-deflection analysis more 
rational# 

♦ The experiment work needs to be carried out to 
support findings of the theoretical investigation such as the 
present one. 
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